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ABSTRACT 


In  this  paper,  a  canonical  formalism  has  been  developed  for 
the  description  of  the  negative  mass  instability  (N.M.I.)  and  longi¬ 
tudinal  oscillations  of  relativistic  beams.  This  formalism  has  been 
applied  to  ascertain  the  stabilizing  effect  of  betatron  oscillations,  and 
to  determine  the  dispersion  relation  governing  counterstreaming  ions 
and  relativistic  electrons.  The  results  show  that  only  the  spread  in 
P  ,  the  canonical  angular  momentum  of  the  particles,  contributes  to 
stability.  The  N.M.I.  equation  for  two  streams  is  the  same  as  though 
each  were  separately  present;  and  the  dispersion  relation  for  longi¬ 
tudinal  oscillations  of  beams  in  a  magnetic  field  is  given  by  the  N.M.I 
dispersion  relation,  and  not^  by  the  dispersion  relation  for  longitudinal 
oscillations  of  collinear  beams.  Moreover,  the  dispersion  relation  for 
longitudinal  oscillations  of  thin  collinear  beams  differs  from  the  usual 
equation  by  a  non-trivial  factor. 
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The  equations  in  the  four  chapters  are  labelled  by  number, 
e.g.  ,  eq.  2-34,  the  34th  equation  in  Chapter  II. 

The  equations  in  the  Appendices  are  labelled  by  letter,  e.  g.  , 
eq.  D-12,  the  12th  equation  in  Appendix  IV. 

Within  the  chapter,  eq.  2-34  is  referred  to  simply  as  eq.  34. 

Within  the  appendix  ,  eq.  D-12  is  referred  to  simply  as  eq.  12. 

The  symbols  used  and  some  words  are  explained  and  defined 
in  the  Definitions  and  Symbols  Section  (p.123  ) 
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CHAPTER  I  -  INTRODUCTION 
Sec.  1. 1  -  The  Description  of  the  N.M.I. 

The  N.M.I.  causes  azimuthal  clumping  of  intense  beams  in  any 
device  with  a  magnetic  field  of  axial  symmetry  that  provides  radial  and 
axial  focusing,  e.g.,  synchrotron,  betatron,  and  mirror  machine 
(e.g.  ,  DCX).  (See  Fig.  1  . )  In  the  unstable  regime  this  clumping  grows 
so  that  eventually  transverse  space  charge  forces  exceed  the  focusing 
forces  and  part  of  the  beam  is  lost.  (A  crude  analogy  is  the  loss  of 
water  in  a  shallow  circular  trough,  when  waves  are  set  up. )  In  the  stable 
regime,  the  clumps  travel  with  the  beam  and  move  very  slowly  relative 
to  it.  Thus  l  equally  spaced  clumps  in  a  beam  moving  at  cyclotron 

frequency  0.  ,  will  give  rise  to  a  charge  variation  at  one  point  of  fre¬ 
quency  U/  *  2  6o  . 

A  simple  explanation  of  the  cause  of  the  instability  is  the  following. 
Consider  an  azimuthally  uniform  distribution  of  cold  particles  in  a  beam. 
Making  a  small  sinusoidal  perturbation  in  the  beam  density  effects  a 
sinusoidal  electric  potential,  rotating  with  the  beam.  Those  particles 
ahead  of  the  potential  bump  will  be  speeded  up  and  those  behind  it  will 
be  slowed  down.  Thus  one  would  expect  the  bump  to  evanesce.  However, 
those  particles  which  were  speeded  up  move  outward  radially  due  to  the 
centrifugal  force.  The  amount  of  radial  motion  depends  on  the  magnetic 
field  shape.  For  weak  focusing  machines  (o  <  h  <  /  ),  the  radial  ex¬ 
cursion  is  large  enough  to  overcompensate  the  increase  in  linear  veloc¬ 
ity  so  that  the  angular  velocity  decreases.  On  the  other  hand,  those 
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particles  that  are  slowed  down  in  linear  velocity,  move  inward  radially 
and  thus  speed  up  in  angular  velocity  so  that  the  net  effect  is  that  the 
particles  tend  to  move  toward  the  angular  position  of  the  potential  bump 
and  the  perturbation  grows.  Since  the  angular  acceleration  is  opposite 

to  the  force,  the  effect  is  as  though  the  particles  have  negative  effect¬ 
ive  mass. 

As  the  radial  focusing  increases  in  strength,  the  radial  ex¬ 
cursion  becomes  relatively  smaller  so  that  an  increase  in  linear  veloc¬ 
ity  causes  an  increase  in  anguLar  velocity  and  the  effective  mass  be¬ 
comes  positive.  This  is  the  situation,  for  example,  in  strong  focusing 
machines  (below  transition  energy)  where  the  N.M.I.  will  not  occur. 

Gravitational  forces  are  weak,  in  the  above  sense,  so  that 
particles  have  an  effective  negative  mass.  Such  a  situation  exists  in 
Saturn's  rings.  There,  however,  the  forces  between  two  particles  are 
attractive  so  that  the  negative  mass  prevents  clumping  and  the  ring 
system  is  stable.  This  was  first  pointed  out  by  Maxwell/  ^ 

The  cause  of  the  instability  may  also  be  seen  by  examining 
normal  synchrotron  operation.  During  the  acceleration  cycle  of  a  syn¬ 
chrotron,  an  R.  F.  field  is  applied  across  a  gap.  If  we  assume  that 
the  gap  is  so  small  that  the  time  change  of  the  gap  field  is  negligible 
while  it  is  being  traversed  by  the  particle,  then  the  energy  gained  by 
the  particle  is  dependent  only  on  its  phase  relative  to  the  R.  F.  field. 

In  this  case  if  the  R.  F.  peak  gap  potential  is  e  V  ,  then  the  potential 
may  be  replaced  by  an  equivalent,  continuous  rotating  potential  over 
the  whole  path  of  the  particle  given  by  JiX  0<X-  ^6  -  ^  ft)  ^ 


The 
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energy  gained  in  one  cycle  of  the  particle  is  equal  in  both  cases. 
Analysi  ■»  of  the  phase  motion  shows  that  there  is  a  point  of  phase 
stability  near  the  peak  of  the  rotating  field.  The  beam  particle  den¬ 
sity  is  therefore  a  maximum  at  the  same  point.  If  the  magnetic  field 
is  stationary  then  the  stable  point  is  at_  the  peak.  Suppose  now  that 
there  is  no  rotating  external  electric  field  but  a  rotating  internal 
field  caused  by  6  perturbations  in  the  beam  density.  Since  the 
stable  phase  point  is  at  the  field  maximum,  the  particle  density  at  the 
maximum  will  grow,  which  will  make  the  field  still  stronger,  resulting 
in  beam  clumping. 

Sec.  1.  2  -  Background 

This  work  was  motivated  by  the  conjecture  that  the  N.M.I.  is 
the  effect  which  most  severely  limits  the  maximum  currents  allowed 
in  a  plasma  betatron.  For  a  plasma  at  an  initial  temperature  of  3  e.v. 
the  linearized  N.M.I.  theory  predicts  stability  at  a  neutral  beam  den¬ 
sity  where  the  current  is  only  £  ^  amp  of  relativistic  electrons  (see 
eq.  2-47). 

The  plasma  betatron  is  a  device  which  accelerates  a  neutralized 
beam  of  positive  ions  and  electrons  so  that  the  space  charge  limit¬ 
ations  of  ordinary  machines  do  not  apply.  Examination  of  the  equi¬ 
librium  conditions,^  by  including  the  effect  of  the  self-magnetic  field 

as  done  by  Schmidt,  gives  the  limit  <  £„  ,  which  permits  1000 
amperes. 

(4) 

Instabilities  were  first  discussed  by  Budker,  the  originator  of 
the  scheme  of  the  acceleration  of  a  neutralized  beam.  Two  of  these 
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instabilities,  the  two-stream  longitudinal  and  the  transverse  (sinuous) 

,  (5) 

are  also  discussed  by  Finkelstein  and  Sturrock,  who  find  stability 

criteria  much  less  restrictive  than  required  for  the  N.M.I. 

Another  instability  which  might  severely  limit  the  maximum 

,  (6) 

beam  current  was  noted  by  Rosenbluth.  It  arises  when  a  beam  of 
particles  passes  through  a  background  resistive  plasma.  Since  in 
principle,  this  effect  may  be  eliminated  in  a  plasma  betatron  by  care¬ 
ful  design,  we  will  not  consider  it  further 
(7) 

Harrison  has  pointed  out  that  the  two-stream  longitudinal 

instability  severely  restricts  the  maximum  currents  in  a  non-relativ- 

lstic  electron  beam.  As  already  pointed  out  in  a  paper  by  Finkelstein 

(5) 

and  Sturrock,  hereafter  to  be  referred  to  as  F-S,for  relativistic 

beams,  V  <■  aooo  for  stability,  however  for  slower  beams  as  shown 

by  Harrison  the  stability  requirement  is  ]L  >-||  Mq  ,  for  cold  electron 

and  ion  beams.  (  V.  is  the  lineal  stream  density  multiplied  by  the 

classical  electron  radius;  V.  is  the  electron  stream  velocity,  the  ion 

velocities  being  small;  g  is  a  logarithmic  geometrical  factor  of  order 

unity.)  The  theory  of  the  N.M.I.  shows  that  even  if  this  inequality  is 

satisfied  each  beam  must  be  hot  enough  so  that  ~  >  liL9  £!£  >  (v£g" 

c  N  ot  »  c 

(Since  for  usual  t)  values,  and  g  is  a  log  term,  this 

inequality  is  independent  of  the  device  considered.)  We  see  therefore 
that  the  N.M.I.  ,  discovered  independently  by  Nielsen,  Sessler  and 
Symon,  and  Kolomenskii  and  Lebedev,  restricts  the  maximum  cur¬ 
rents  in  a  plasma  betatron  more  severely  than  the  other  effects. 

Plasma  betatrons  have  been  built  by  Budker  and  Naumov^* 
and  by  workers  at  CERN  with  the  result  that  maximum  currents 
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were  10  amperes,  much  below  the  design  value  of  these  machines. 

The  limitation  is  possibly  due  to  the  N.M.I.  A  plasma  betatron  is 

(12) 

also  under  study  at  Stevens  Institute  of  Technology. 

This  instability  is  of  broader  interest,  because,  as  we  shall 

show,  thin  beams  in  mirror  magnetic  fields  where  <><  h  c  I  ,  will  also 

(13) 

be  subject  to  it,  e.g.  in  the  DCX  machine,  where  it  may  be  the 

(14) 

cause  of  the  observed  frequencies,  as  also  pointed  out  by  Fowler. 

(IS) 

The  Astron,  containing  a  beam  in  a  mirror  field  may  also  be  sub- 

(!6)  (17) 

ject  to  the  N.  M.  I.  Samoilov  and  Seidl  have  observed  particle 
bunching  in  betatrons  and  attribute  this  to  the  N.M.I.  However,  suf¬ 
ficiently  detailed  measurements  have  not  been  made  to  verify  this 
conjecture.  These  authors  also  suggest  that  the  N.M.I.  is  the  major 
cause  of  capture  of  particles  into  stable  orbits  in  betatrons. 

Sec.  1.  3  -  Outline 


We  have  derived  the  N.M.I.  equations  using  a  canonical 

formalism  and  the  relativistic  Hamiltonian.  This  procedure  allows 

one  to  include  additional  effects  easily.  With  the  resultant  dispersion 

relations,  we  derive  a  necessary  and  sufficient  criterion  for  stability  , 

which  is  simple  only  for  single  humped  distributions.  The  stability 

criterion  may  be  given  explicitly  for  Maxwellian  distributions,  and  is 

similar  to  the  result  obtained  for  rectangular  pulse  distribution 
,  .  (8),  (9) 

functions  by  other  authors  and  hence  justifies  the  use  of  pulse 

functions.  Our  dispersion  relation  differs  somewhat  from  earlier  re¬ 
sults  so  that  stable  distributions  exhibit  damped  oscillations.  This 
effect  is  shown  explicitly  for  a  resonance  distribution  function. 
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(9) 

Kolomenskii  and  Lebedev  have  obtained  similar  stability  criteria 

for  a  resonance  function  but  have  not  worked  out  the  damped  situation. 

(13) 

This  damping  is  mathematically  analogous  to  Landau  damping  in 
infinite  plasmas.  Our  results  exhibit  no  damping  for  pulse  distribution 
functions  in  agreement  with  the  results  of  other  authors. 

Previous  works  of  other  authors  have  dealt  with  circulating 
beams  enclosed  within  conducting  boundaries.  This  paper  considers 
unshielded  beams  so  that  every  part  of  the  beam  sees  every  other  part. 
We  find  that  for  small  wave  numbers  of  the  perturbation,  in  the  rela¬ 
tivistic  domain,  the  beam  will  be  stable  even  if  it  is  cold.  (If  —  = 

J _ .  *1* 

k*  c»  18  negative  then  there  is  stability,  where  3... 

and g  x  J"  for  typical  cases. ) 


Next  we  consider  the  effect  of  betatron  oscillations  on  the  N.M.I. 
and  treat  separately  the  axial  (z  )  and  radial  (r)  oscillations.  We  find 
that  these  have  a  very  slight  effect  on  stability  and  therefore  only  the 
spread  of  contributes  to  the  stability.  Our  calculations  also 

show  that  the  growth  rate  of  the  instability  slows  down  as  it  approaches 
the  radial  betatron  oscillation  frequency  where  the  equations  break  down. 

Finally  we  generalize  our  equations  to  find  the  dispersion  re¬ 
lation  under  the  N.M.I.  for  two  streams,  counterstreaming  ions  and 
relativistic  electrons.  We  find  that  the  stability  criteria  are  almost 

the  same  as  though  each  beam  were  present  by  itself.  The  difference 

is  that  when  h  >  4-  and  /  >>  3  >  _L  ,  , 

a  j  f  ,  the  electron  modes  are 

stable,  even  for  a  cold  electron  beam.  The  ion  modes,  however,  are 

still  unstable. 
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Next  it  is  shown  that  the  equation  for  longitudinal  oscillations 
of  collinear  streams  is  valid  for  circular  streams,  only  when  « 

-L  ,  which  is  not  generally  true.  (It  is  true  only  in  strong-focusing 
machines  below  transition. )  The  equation  for  longitudinal  oscillations 
may  therefore  be  derived  from  the  two-stream  N.M.I.  equations  by 
letting  — >  O  .  This  equation,  valid  only  for  small  V  because 

the  N.M.I  equations  are  restricted  to  small  V  ,  differs  from  the 

(5) 

F-S  longitudinal  equation.  Going  back  to  the  basic  equations,  an 
equation  for  the  longitudinal  oscillations,  valid  for  all  V  ,  is  ob¬ 
tained  which  gives  stability  for  even  higher  currents  than  found  in  F-S. 
To  check  the  validity  of  our  equation,  the  dispersion  relation  for  two 
infinitely  wide  beams  is  obtained  from  it.  If  the  ion  beam  is  stationary 
and  cold,  the  dispersion  relation  agrees  with  that  found  by  Bludman 

(19) 


et  al.  If  we  set  the  number  of  ions  equal  to  zero,  we  find  that  the 
resulting  dispersion  relations  may  be  obtained  by  a  Lorentz  transform¬ 
ation  from  the  dispersion  relations  of  both  thin  and  infinite  beams  of 

non- relativistic  electrons.  These  results  show  that  our  modification 

/  n  \9- 

of  the  F-S  equation,  which  consists  of  a  factor 


'■(tt)'1 


is  cor¬ 


rect.  This  equation  is  similar  to  cq.  32  (eq.  9  in  the  abridged  trans- 
(4) 

lation)  of  Budker,  who  lias  obtained  the  same  stability  criterion.  (His 
equations  neglect  beam  temperatures  and  the  ion  beam  velocity.) 
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CHAPTER  II  -  ONE  STREAM 


Here  we  develop  the  dispersion  relations  for  a  single  stream. 

A  constraint  equation  is  found  which  reduces  the  Boltzman  Equation  to 
a  one -dimensional  equation.  Since  the  equation  for  the  potential  is  then 
given  by  an  integral  over  one  momentum  variable,  the  resulting  equa¬ 
tions  are  formally  similar  to  the  one-dimensional  system  first  studied 
(18)  (20) 

by  L.  D.  Landau  and  later  also  by  Backus.  These  equations  are 

(21) 

solved  following  Jackson  while  the  Nyquist  diagram  technique  of 

(22) 

Penrose  is  used  to  obtain  generalized  stability  criteria  including 
the  stability  criterion  for  a  Maxwellian  distribution.  These  equations 
are  then  solved  exactly  for  a  resonance  function  and  a  pulse  function. 
Finally  the  Boltzman  Equation  is  solved,  non-relativistically ,  by  in¬ 
cluding  the  (z)  axial  betatron  oscillations  and  the  (r)  radial  betatron  os¬ 
cillations  separately. 

Sec.  2.1  -  The  Constraint  Equation 

The  basic  equation  for  our  system  is  the  collisionless  Boltzman 
equation 

if  V  e  If  -  fa  ,  r  if .  fr  Of,  j,  Df  2-1 


d9 


t)  /Jo 


dr 


2/’, 


(9 


If  the  coefficients  cj.  and  p,  are  obtained  from  a  relativistic  Hamiltonian 

and  describes  particles  with  the  same  rest  mass,  then  this  equation 

is  relativistically  correct.  A  brief  discussion  of  the  relativistic  invariance 

(23) 

is  given  by  Belyaev  and  Budker.  In  the  following  way  they  show  that 
^  is  a  Lorentz  invariant  scalar.  The  particle  flux  and  density  four- 
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v  cttur, 


jk  is  obtained  from  J  Fk  d ^ 


,  invariant  expre* 


siom  In  terms  o£  the  four-velocity  U„  ,  F„  =  FW(,  and  F  (V.  ,  fh  )  = 

^  •  ft-  *  ^  (  H?"  ’  |  ^t)J  ~b..cy  The  5  function  arises  because 

F  describes  particles  with  the  same  rest  mass,  f  i.  the  attribution 
function  of  eq.  1,  Bo  that  Jfj>f  „  „  ,  the  number  of  particle8  per 

unit  volume.  Since  the  ,5 function  ia  written  in  an  invariant  way.  ¥~ 
muat  be  invariant  becauae  F  is  invariant.  (There  is  an  error  in  the  ex¬ 
pression  for  F  and  H  in  their  paper. '  The  factor  in  F  should  be  -i  , 
not  ic  ,  and  H  should  be  multiplied  by  c  .)  Because  of  the  S  function. 
¥  a  function  of  only  seven  variables,  the  four  coordinates  and  the 
three  momenta.  Eq.  1  may  be  derived  from  their  invariant  Boltaman 
equation  by  integrating  i,  over  J  f,  to  eliminate  the  J  function.  Thus 
eq.  1  is  relativistically  correct. 

Our  procedure  is  to  simplify  the  expressions  for  the  coefficients 

of  eq.  1  and  then  solve  the  Boltzman  equation  by  a  perturbation  procedure. 
This  means  writing  f  =  £  ♦  jf  ,  where  £  «  £  in  some  operational 

sense  and  £  describes  the  unperturbed  configuration  which  is  time  in¬ 
dependent.  If  ,  initially  small,  has  an  exponentially  increasing 

time  dependence  then  the  system  is  unstable. 

The  unperturbed  system  consists  of  one  specie  of  particles  ro¬ 
tating  about  an  axially  symmetric  magnetic  field.  (See  Fig.  i  ).  The 
Bz  field  falls  off  slowly  with  radius  near  rQ  according  to  Bz  =  B0  )" 
where  Oi  h  el  to  provide  focusing,  as  explained  in  Appendix  I.  The 

particles  occupy  a  toroidal  region  of  small  cross  section  and  form  an 
azimuthally  uniform  distribution. 


i. 
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The  equations  for  the  transverse  motion  are  derived  in 
Appendix  I  and  are  given  by  eqs.  A-8,  9  and  10,  which  describe  the 
usual  betatron  oscillations.  Note  that  m,  the  relativistic  mass 
as  defined  in  Appendix  I)  is  constant,  since  the  energies  of  betatron 
oscillation  are  constant.  These  equations  for  the  transverse  motion 
are  valid  if 


bK  (  dti 

3r  Sj- 

are  neglected  relative  to 


3  <p°  3  <p' 


hr 


hi 


2-2 


in  eqs.  A-3a,  3c.  (The 


h_A_a  d  A  $ 
dr  '  ht 

superscript  zero  refers  to  quantities  due  to  the  unperturbed  beam.) 

I  -  h  /l  '  j  v a 

T~ 


This  requires  that 


or  ion^  stream.  It  is  also  required  that  Ar,  Az  and  A  be  negligible. 

In  fact  Ar  =  Az  =  0  because  of  the  symmetry  of  the  particle  motion. 

(3)  W 

Finally,  as  Schmidt  has  pointed  out,  the  self-field  term, 

a  r 

causes  a  radial  shift  in  equilibrium  orbit  which  is  negligible  if  g 

in  which  case  Au0  is  also  negligible. 

The  Schmidt  criterion  may  be  derived  in  the  following  manner. 

0 

Bz  =  0  at  the  center  of  the  current  torus,  i.  e.  ,  near  r  =  rQ,  z  =  0, 

if  the  current  is  distributed  uniformly  over  the  cross  section.  The  re¬ 
lation  Bz  =  ^ 6  *  Bp  -  Q 


M*'*- 


for  the  electron 


(-) 


then  gives 


F  rom 


__  hAl  ...  Ae 

-dr  r  W  ~ 

Appendix  I  we  see  that  the  radial  motion  of  the  particles  is  determined 

by  a  vector  potential  Ag  .  The  beam  center  is  at  the  bottom  of  the  well 

defined  by  the  total  vector  potential,  the  external  plus  the  self-field. 

See  Fig.  2.  The  location  of  the  bottom  is  given  by  the  solution  of 

=  O  Or  lA  ■  As  using  Bz  =  0.  To  find 
dr  hr  ar  <J r  r 
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the  value  of  Ae  at  the  center  of  the  beam,  near  r  =  rQ,  we  use  eqs.  A-5b, 
C-17  and  C-20  to  obtain  K  “  ll’  IL  ^  -  V'jo  [Te  .  By  eq.  A-2b, 

:  -  iAl  --  (A* +  /a) 

>H,JC 

solving  for  A°d  gives  A°e  ' 

(1  .  n)  klLt.  At 


Mt 


J  to 


bottom  of  the  well 

( l-h )  LlL*.  Aai 

rS 

reduces  to 


.3  TT  C 

.  Substituting  this  value  of  \fg  and 

"±A  he 

r  -r  .  By  eq.  A -4, 

1 — +  ‘"I*  3r 

.  -y  — 

•  Thus  using  the  equation  for  the  location  of  the 

r  jr  ,  we  obtain  finally  —  — h  —  — = 
d  '  yri+^% 

.  Since  )L3j  «  I  ,  and  r  *  rD,  this  expression 


The  maximum  of  the  R.H.S.  of  this  equation  is  ,  which  gives  us 
Schmidt's  criterion.  Note  that  the  current  loop  moves  radially  inward 
in  the  Betatron  field,  contrary  to  a  free  current  loop  which,  as  is  well- 
known,  expands. 

The  unperturbed,  zero-order  azimuthally  symmetric  distribution 
describes  particles  with  a  spread  in  pe  values  and  a  range  of  betatron 
oscillation  amplitudes. 

The  variables  p^  ,  r  -  rQ,  and  z  are  considered  first-order 
small.  Quadratic  terms  in  these  quantities  will  be  neglected.  As  a  re¬ 
sult  of  the  perturbation  which  causes  azimuthal  fields,  pc  is  no  longer 
constant  for  each  particle  but  changes  slowly  with  time  (see  eq.  A-3b) 
and,  therefore,  m  will  too.  There  will  now  be  terms  due  to  Ar,  Ae 
and  <p  ,  in  the  expression  for  pz,  and  hence  additional  terms  in  eq. 

A-9 .  Az  remains  zero  because  the  motion  is  symmetric  about 

the  z  -  0  plane.  Ar  may  also  be  neglected  for  thin  beams  because  the 


radial  phase  velocity  due  to  the  N.M.I.  is  always  much  slower  than  c. 
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The  terms  3/0  ,  ,  3^?',  give  the  effect  of  the  perturbed 

&r  ’  3?-  3 r  7?" 


transverse  space  charge.  It  appears  plausible  that  under  the  conditions 
that  the  zero  order  transverse  space  charge  effects  may  be  neglected 
relative  to  the  focusing  betatron  field,  that  the  perturbed  terms  may 
also  be  neglected.  It  is  possible  in  fact  to  show  this  non-relativistically 
with  the  formalism  of  Sec.  2.  5.  This  suggests  that  if  V£|  a  these 
terms  may  be  neglected. 

With  the  above  assumptions,  we  obtain  from  eqs.  A-8  and  A-9 
that  pz  =  mi  and  pz  =  -  b  i  .  Since  m  must  be  constant  for  particle 


motion  in  a  static  magnetic  field,  these  equations  give  mz  =  -  h  i  or 
pz2  +  bz2  =  constant.  This  means  that  the  energy  of  the  z  betatron  oscil¬ 
lations  is  constant.  Thus  if  is  a  function  of  z  and  pz  only  through 
pz2  +■  bz2,  i.e.  ,  Y  {r,  p  ,  0  »  P,  .  t,  pz2  +  bz2),  the  two  z  terms 

in  the  Boltzman  equation  add  to  zero  as  may  be  verified  by  substitution. 
The  coefficients  of  the  other  terms  do  not  contain  z  or  pz  to  ^rst  order 
so  that  we  may  integrate  the  Boltzman  equation  over  c >  and  writing 
Y~  f  Yd  id  fr  obtain, 


6  *  foil-'.  r  2-3 

d  t  3  0  3  fo  d  fr  3  r 


We  have  thereby  reduced  the  equation  to  a  two-dimensional  one. 

Next  we  deduce  a  constraint  equation  linking  r  and  pff  ,  through 
which  the  problem  is  reduced  to  only  one  dimension.  The  equation  of 
motion  in  the  r  direction  is,  by  eq.  A-10  (again  neglecting  the  transverse 

t 

space  charge  forces  and  Ar), 

d_  wi  f  -  -  (-  b; )  3_  f-  fp  t  A  a ) 

dt  VI  r  c  >  2_4a 
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For  m  and  p4  constant  this  equation  may  be  rewritten  as 

-A  ft  -  K*(r-n) 


2-4b 


r0 


where  ior  -  ji-h  =  JT7 77 ^ - e  fi.j  .  This  shows  that  a  particle  oscil¬ 
lates  about  an  equilibrium  position  r*  =  r„  +  —  _£ _  fJn  with 

a  frequency  to/r  .  This  equilibrium  position  will  vary  with  p6  ,  but 
the  frequency  remains  constant  to  lowest  order.  This  may  be  seen 

more  clearly  by  examining  the  plot  of  Afi  in  Fig.  2.  Eq.  4a  gives  the 
motion  of  a  particle  in  the  potential  well  I?  -  £  A  B  .  Clearly  the 
minimum  of  the  total  well  is  shifted  according  to  the  value  of  p&  . 

(This  is  seen  by  adding  the  curve  (--£■  f0  )±  to  AB  in  Fig.  2). 

Under  the  influence  of  azimuthally  varying  electric  fields,  p  will 
change,  as  appears  from  eq.  A-2b.  These  fields  will  occur  as  a  result 
of  the  N.M.I.  Suppose  now  that  the  p0  variation  of  a  particle  is  very 
slow.  Then  if  the  particle  is  initially  at  the  bottom  of  the  total  well,  it 
will  stay  very  near  the  bottom  and  follow  the  shifts  in  p0  .  This  may  be 
shown  easily  by  writing  the  steady  state  solution  of  eq.  4  with  the  initial 
condition  that  the  particle  is  resting  at  the  bottom  of  the  well,  i.  e.  , 

r=r=0  at  t  =  0  and  the  assumption  that  p&=  pcoa  .  This 

solution  is 


r-r„ 


-  6t 


^ 


U‘ 


2-5 


^  To  ^w3  -  Ur3 

If  now  the  pfl  oscillation  is  so  slow  that  ^  ,  then  the  second 

term  on  the  R.  H.S.  of  eq.  5  may  be  neglected  and  we  find 


r  -  r„ 


/-•  h 


*  So  r. 


Po 


2-6 
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i.e.  ,  the  particle  follows  the  bottom  of  the  well.  This  may  also  be 
seen  from  eq.  4b  as  r  is  now  much  less  than  either  term  on  the  R.H.S. 
of  the  equation.  Setting  r  =  0  gives  the  same  eq.  6.  This  is  the  de¬ 
sired  constraint  equation.  This  means  that  if  the  particle  is  initially 
at  the  bottom  of  the  well,  with  nearly  zero  amplitude  betatron  oscil¬ 
lations,  then  if  the  bottom  of  the  well  shifts  slowly  enough,  the  particle 
will  follow  the  bottom  without  any  fast  betatron  oscillations  being  ex¬ 
cited.  Their  amplitude  by  eq.  5  is  only  the  amplitude  of  the 

u/  1 
r 

slow  motion. 

Due  to  the  p0  changes  and  the  azimuthal  field,  the  relativistic 
mass  m  will  change  with  time,  so  that  strictly  a  term  should 

be  included  on  the  L.H.S.  of  eq.  3.  This  term  is  however  of  order  p* 
and  is  neglected  as  we  keep  only  terms  -^p#. 

The  restriction  uj  a.  k/f  ,  naturally  places  restrictions  on 
the  solution  of  the  Boltzman  equation  describing  the  N.M.I.  which 
restrictions  we  now  derive.  We  shall  assume  below  that  ‘ 

The  particles  travel  at  an  average  velocity  0o ,  so  that  6=6et.  Hence, 
the  time  variation  of  p#  for  a  particle,  is  f0  ~  (  and  the 

frequency  of  oscillation  of  p9  of  a  particle,  is  2©.  -Jl  .  Thus  the 
above  condition,  that  the  constraint  equation  be  valid  is 


«  I 


or 


Si  -  U>. 


I  l-b  &• 


«  I 


1-1 


which  may  be  verified  to  be  consistent  with  the  dispersion  relation 
obtained  below. 
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We  shall  now  show  how  the  constraint  equation  may  be  used  to 
reduce  the  Boltzman  equation  by  one  more  dimension.  Consider  the 
quantity 

x  2  r  -  n.  =  r  -  rD  *  —  — -  fo 

i  -  *  e  ro 


which  measures  the  deviation  of  the  particle  from  the  bottom  of  the 
well.  We  shall  make  a  transformation  of  from  the  variables 

6,  t,  ,  pr,  r  to  Q  ,  t,  pffI  pr,  x.  Keeping  in  mind  the  fact  that 
X=  -f(ri  ft )  we  may  write  the  last  three  terms  of  the  Boltzman  equation, 
eq.  3,  as 


d  fa  3  fr  9  r  \2  fo  c)  f6  3  x  )  2  fr  3  X 


-Or-'  -\X^I 

where  I*  ~  1*  { 6 ,  t,  pfl ,  pr  ,  x).  We  may  now  substitute  for  the  co¬ 
efficients  p  and  r.  Neglecting  again  K  ,  we  obtain  r  -  from 

r  hi 

eq*  A- 2a  and  f>r  =  -tnu/fx  from  eq.  4a.  Thus  we  get  for  the  last  two 
terms  of  the  Boltzman  equation 


rl£'- 

3  fr  2  X  3  /V 


2-8 


« t;  -  r  t,  p0 ,  pr^  +  m^  w2  x2),  eq.  8  equals  zero  so  that  the 
two  r  terms  now  give  zero  in  the  Boltzman  equation  which  then  be¬ 


comes 
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IT-'  *  e  iT'  *  /ii  lT')=  o 

3  30  la/3©  3  3  x  / 

We  shall  now  integrate  this  equation  over  dn  d pr  .  Thus  the  first 

term  gives  £  Z  ,  where  In  the  next  term 

8  f- 

we  must  be  more  careful  as  0  =  .  We  will  assume  that  the 

functional  dependence  of  ^  on  is  sharply  peaked 

about  x*  =  0  .  The  derivation  of  the  constraint  equa¬ 

tion  shows  that  if  pf  and  x  are  zero  initially  that  they  remain  very 
small.  Thus  X  can  be  a  sharply  peaked  function  of  these  vari¬ 
ables.  Setting  x  =  0,  now  means  that  6  ~^(fo)  only  because  the 
constraint  equation  is  valid  and  r  is  a  function  of  p#  .  Thus  integrating 
the  second  term  we  obtain 


6 


dir" 
*  e 


where  6-  f(fe )  now.  The  fourth  term  is  odd  in  x,  because 

3/  3  r;  3X  .  dt/ _ •  2 

3  fd  3  x  3  fo  3  (/>,*  *  *»’  •*?*)  3  x 

The  first  factor  is  a  constant.  The  second  is  even  in  x,  while  the  third 
is  odd.  Thus  the  integral  over  dx  gives  zero.  This  leaves  only 

3  ir‘ 

fo  _  _  - 

2  fo 

Integrating  this  term  now  over  dx  dpr  now  gives 


Til 

_ 

3  fo 

where  again  we  must  use  the  constraint  equation  to  eliminate  any  r 
dependence  in  p^  . 
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The  validity  of  the  reduction  of  the  three-dimensional  Boltzman 
equation  to  a  one-dimensional  equation  is  justified  more  rigorously  in 
Sec.  2.5,  where  the  complete  solution  of  the  Boltzman  equation  is  per¬ 
formed. 

If  the  magnetic  field  does  not  satisfy  the  Betatron  2-1  condition, 
but  still  satisfies  0  *  h  *  /  ,  then  all  the  results  presented  in  this 

paper  are  still  correct,  because  as  shown  in  Appendix  II,  in  such  a 
system 

fo  _  £  /I  *  -  f±  -  !L  A  b 
r  r 

where  represents  the  external  field  (i.  e.  ,  mirror  fields  or  synchro¬ 
tron  fields)  and  Pj  =  f>0  -  f0j  .  (f0o  corresponds  to  an  equilibrium  orbit 
at  rQ. )  Therefore  if  Pj  replaces  p#  in  all  the  equations  of  this  paper, 
they  will  still  be  correct  because  all  the  equations  of  Appendix  I  are 
the  same. 

Sec.  2.  2  -  The  Basic  Equations 

The  Boltzman  equation  for  the  system  is  now  one-dimensional. 
Setting  =  f'  ,  we  have 

+  e  li'  i  f>0  ll.  =  o  2-9 

dt  d  6  3  fo 

where  the  number  of  particles  in  an  element  ciodfo  is  given  by 

d  N  -  2-10 

while  the  coefficients  6  and  p6  are  defined  by  eqs.  A -2b  and  A- 3b, 

-  .  l  t±  -  LAo 

*  h-i^r  \  r  C 


e 


2-11 
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■  _  e_r  34r  4  <2  4  3^  +  ej£<  -  e 

c  29  c  2e  c  J $  de  2-12 

Eqs.  9,  11  and  12  are  still  too  complicated  to  be  solved  exactly.  We 
shall,  instead,  use  a  perturbation  expansion,  V/=  V#  *  ^  where 

describes  the  time  independent,  azimuthally  uniform  distribution, 
while  ^  ,  contains  the  d  and  t  dependence,  and  is  a  small  quant¬ 
ity  compared  to  'H,  .  This  is  consistent  with  our  previous  approxi¬ 
mations.  Thus  O  ,  while  -f0  ~  $  »  by  eq.  12. 

Thus  to  terms  of  lowest  order,  the  Boltzman  equation  now 
becomes 

Hh  f  9  +  f*  “  O  2-13 

S  t  2®  3  fo 


This  is  called  the  linearized  equation  because  all  terms  are  linear  in 
,  .  ^  d> 

r,  .  Note  that  f„  — — -  has  been  dropped  as  it  is  of  second 

o  fo 

order  in  ^  .  The  two  first  terms  in  eq.  12  may  also  be  dropped 

as  they  are  of  higher  order  than  the  last  two.  More  particularly, 

Az  =  0  because  the  motion  in  the  z  direction  is  symmetric  so  that 
Iz  =  0.  The  first  term  may  be  neglected  because  by  the  constraint 
equation,  eq.  6,  r  ~  p&  ,  while  Ar  ~  ,  so  that  this  term  is  second 

order  in  ^  .  The  coefficients  Q  and  p^  may  now  be  written  more 

explicitly.  They  are 


y  bur 


2-14 


22 


where 


and  *>a  «  2-15 

'  c  2e  dO 

Since  9  multiplies  ^  ,  in  eq.  13,  a'„  gives  a  second  order 

term  and  is  neglected.  The  rest  of  0  is  a  given  function  of  p0  and 
r,  and  through  the  constraint  equation,  T-  f(fo)  ,  is  a  function  of  ps 
only.  (Ag  is  defined  by  eq.  A-4  where  now  z  -  0.)  Since  pa  is  small, 
6  may  be  expanded  as  a  linear  function  of  pfi  as  detailed  in 
Appendix  IV.  The  result  is 


6  *  6.  -  k  f0 


k 


_ ! _  f—L 

*-  r 1  i  -  * 


j_\  2-16 


This  gives  one  coefficient  of  the  linearized  Boltzman  equation, 
eq.  13.  The  other  coefficient  is  p9  .  We  desire  its  explicit  depend¬ 
ence  on  pe  and  V?  also.  It  now  proves  more  convenient  to  use  Fourier 
and  Laplace  transforms  as  defined  by  the  expressions 


C  =  f  Ji  (  die  e 

\  K  l 


fo 


fo(o+) 


2-17 


which  imply  the  reciprocal  relations 


toW)\ 


I 

(ait  j1 


^  i  lo 
£  e 


i-  -<* 


inf 


2-18 


23 


where  the  contour  W  is  chosen  in  the  upper  half-plane  parallel  to  the 

real  Si  -axis  above  any  poles  in  ^?/l  »  or  Pe  •  oee  FiS-  3  . 

To  find  p  we  consider  first  c/>  ,  whose  transform  is  given  as  a 
ro 

function  of  ^  by  eq.  C-20 


V5  ?/l  =  e  3i  ft  dfo 


and  is  valid  only  in  the  beam,  and  when  the  wave  length  of  the  azimuthal 

perturbation  is  much  larger  than  the  beam  width.  Since 

and  'A  gives  rise  to  an  azimuthally  symmetric  f  -  djj 

do  3© 

and  by  eq.  19 


iur  ■  *r 


Next  we  examine  the  expression  for  Aq  .  F rom  eq.  A-5 

^  _  r  \Te  ooa(fl-g’)  j>  (T,o)  Sr  r  dp 

6  J  ~  Ir-r'l 

2-21 

_  e  /  Vo  Cou-  (e  -  Q ')  ^(0,  fo )  do  dfo 

)  c  jr  - r'l 


because  the  component  of  I  along  Aa  is  I  cos  ( 0-9  ).  This  expression 

for  A q  will  be  substituted  into  eq.  15  for  p0  and  hence  must  be  evalu¬ 
ated  to  first  order  in  */?  .  Higher  order  terms  will  be  dropped. 

We  need  first  an  explicit  expression  for  v'e  =  / (fo)  »  valid  to  first 
order  in  ,  where  V~0  ~  r  Q  •  For  r,  the  constraint  eq.  6  gives, 


r  -  r 


fo 


2-22 


24 


the 


In  the  expansion  of  6  in  Appendix  IV,  some  changes  must  be  made 
as  we  desire  9  to  first  order.  We  now  include  Ae  ,  so  that  the 
expansion  will  have  ll  -  as  a  factor.  Also  we  will  expand 

r  dependence  separately  to  lowest  order  in  r  -  rQ  and  then  substitute 
eq.  22.  Thus  using  eq.  14  instead  of  eq.  D-l  and  D-3  in  Appendix  IV 
we  obtain 

•  *o. 

and  to  first  order  terms, 

v.»  ro  =  (r,  fa  +  J  r>yAr*  'r°rA^' 


r  0  +  __L__  /  f°  ~  r°  —Ae  ) 

•e-  kur.'i’l  ‘  ' 


2-23 


=  t/; 


__L  (£t  'LA') 

i-’hi.  (  r.  c  ' 


This  term  together  with  K)/-  Vi  *  ^  >  must  be  inserted  in  eq.  21. 

Keeping  only  terms  first  order  in  7, 1  ,  we  obtain  from  eqs.  21  and 

23 


i  f» _ W  ~  G  A'  ^ 


Coa.  (e  o') 


dedf6 

/  r -r*  / 


2-24 


The  $  dependence  of  the  bracket  is  only  in  and  A1,  as  A1 
The  results  of  Appendix  III,  in  particular  eq.  C-20,  now  gives 


nn  _  e  q, 

0 

To  ) 


(\ 


*  t 

c 


in 


,<  /  n 


h..  c  s-‘rJ 


^ .  c  > 5 


&  fo 


2-25 


25 


The  ^  integration  may  be  easily  performed.  Since  w  is  normalized 

to  N  on  the  field  6  and  pfl  and  is  constant  from  6  =  0  to  2tt,  j J f9 

-  where  N  is  the  total  number  of  particles.  Thus  transposing  and 
<3  TT 

dividing 


>3 

Substituting  20,  23  and  26  into  eq.  15  now  gives, 


As  we  will  show,  the  pe  terms  in  eq.  25  may  be  neglected. 
Eq.  27  is  then  proportional  to 

A 

,  U  ^  1  I  /,  f  v5\  -  Sti  -6  '  iis  “ 

'-(fj  T,,_  j  j  -j-  Pc-,  v 
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(where  J‘ tj  -  -  j  >  0  and  /  from  the  Schmidt  criterion  for 

;  . 

the  neglect  of  the  zero-order  self-field.)  The  term 
the  integral  is  negligible  if  fo  <c  ^in.cr  or  A  t/l  <:< 


outside 


c_ 

J J 


while  for  the  positive  specie  A  v't  «  C  is  sufficient,  if  the  ions  are 
non-relativistic.  Now  consider  the  pa  term  under  the  integral.  By 
eq.  31,  t  consists  of  two  parts.  If  the  part  proportional  to  R  (2). 
(eqs.  30  and  31),  is  sufficiently  smooth  and  involves  only  small  p0  , 
so  that  f0  «  >■  c  ,  (  fu  <-1  3”,  c  r„  for  the  ions)  this  part  may  be 


26 


ne 


glected.  Using  the  other  part  of  ,  and  the  identity  of  eq.  3-32, 

Q  -  JJ- 

may  be  taken  out  of  the  integral  and  equals  - L  .  Thus  sub¬ 

stituting  for  k  from  eq.  16,  if 


6  -  ^ 
'  ? 


«  —  this  term  may  be 

l-h  r„ 


6  -  — 
f  l 


/  -  h  To 


neglected.  For  the  positive  specie,  the  criteria 
is  obtained,  which  is  well  satisfied,  because  the  constraint  equation 
must  be  valid  and  hence  eq.  7  must  be  satisfied. 

Thus  eq.  27  becomes 


•>n  x 

fo  =  ~±JLl 
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±  B  1  - 
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This  equation  together  with  eqs.  13  and  16  are  completely  equivalent  to 

(18) 

the  plasma  system  first  solved  correctly  by  L.  D.  Landau. 

Sec.  2.  3  -  The  Dispersion  Relation 

We  solve  the  equations  using  Fourier  and  Laplace  transforms 

(21) 

following  J.  D.  Jackson's  notation  and  method.  We  multiply  eq.  13 
-3  7T  -'t(lO-Jli) 

by  /  de  I  cit  e  to  obtain  after  integrating  by  parts, 
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On  the  assumption  that  fL  has  a  positive  imaginary  part,  the  first 
term  evaluated  at  t  =  oo  vanishes  and  the  rest  is  -R  (l)  where 

RU)--  f'v  6 

Jo 
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i.  e.  ,  the  Fourier  transform  of  the  initial  displacement.  With  the  defin¬ 
itions  of  the  Fourier  transforms  as  given  by  eqs.  17  and  18,  eq.  29  now 


gives  the  result 


R(l)  -  ^ts.  f*  )  — rrx 

{  '  Z  fo  )  *(2o  -Jt) 


From  eq.  17,  by  integrating  by  parts 


and  thus  eq.  28  becomes 


if. 


Inserting  eq.  31  into  eq.  32,  now  gives 

e’l  (  f>"  -  *0) 


r„  r*  ?  fo 


Since  p^n  is  independent  of  pfi  ,  we  can  solve  eq.  33  for  pe  , 


(f L£I2_  j, 
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where  <])  and  H  are  the  numerator  and  denominator  of  eq.  34  respect- 
ively.  The  inverse  p0  ,  which  gives  the  time  behavior  of  p^  is  ob¬ 


tained  from  eq.  30  and  is,  using  eq.  34 


28 
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Since  the  curve  W  may  be  closed  in  the  lower  half-plane  of  Jl  (see 
Fig.  3  ),  we  may  evaluate  the  integral  using  residues,  and  obtain 


fo  - 
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where  the  sum  is  over  the  poles  of  the  term  in  brackets.  If  the  initial 
perturbation  is  sufficiently  smooth,  then  $>(7,1)),  will  not  contribute 
any  poles  and  the  poles  will  occur  only  for  the  zeros  of  the  denominator. 
If  the  poles  have  a  positive  imaginary  part,  then  f>J  ~  g**  antj 
the  solution  is  unstable.  Since  the  above  function  H  ^r)  is  defined 
only  for  JL  with  positive  imaginary  part  as  appears  from  eq.  29,  we 
must  find  the  analytic  continuation  of  H (y)  in  the  integrand  of  eq.  35 
in  order  to  find  the  residues.  To  find  H  (j-)  explicitly  we  must  insert 
the  value  of  0  from  eq.  16.  We  then  see  that  is  not  continuous 

across  the  real  axis  (viewed  as  a  function  in  the  complex  Jl  plane), 
because 


Z  ( 

~  -j>  U-  1  (  t 


I 
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The  /  function  gives  the  discontinuity.  This  shows  that  H,  defined  by 
eq.  34,  is  discontinuous  across  the  real  axis  with  a  jump  equal  to 

aH  =  H  Ui  U)  -  =  u-i6) 


2  7 ti  e*<]  d  {-Hie. 
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where  we  have  used  the  J  function  properties  S  (/)  =  /(-*),  $(<**)  -  /fcjand 

CL 

recalled  that  the  integration  variable  in  H  is  p&  .  Since  the  analytic 
continuation  of  H  must  be  continuous,  we  add  AH  to  H  defined  by 
eq.  34  to  get  the  form  of  this  function  valid  in  the  lower  half-plane.  The 
dispersion  relations  are  now  given  by  setting  H  =  0.  Thus 

'b  V* 

h$=  LIK 
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In  eq.  39  it  has  been  assumed  that  7>o  .  If  2  *  o  then  it  may  be 
seen  that  the  first  equation  is  unchanged,  but  that  in  the  second  there 
is  a  plus  sign  in  front  of  the  2iri  in  the  third  term  on  the  right  hand 
side. 


Sec.  2.4  -  Stability  Criteria  and  Dispersion  Relation  Solutions 


In  this  section  we  will  summarize  some  results  obtained  by  in¬ 
vestigating  the  solution  of  the  dispersion  relations.  If  one  desires  to 
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know  only  stability  criteria,  i.e.  ,  the  condition  that  It  in  eq.  18  has 
a  negative  imaginary  part,  so  that  the  corresponding  Fourier  compon¬ 
ent  is  exponentially  damped,  then  it  is  unnecessary  to  solve  the  dis¬ 
persion  relations  completely.  By  means  of  the  Nyquist  diagram  as 

(22)  (21) 

elaborated  by  Penrose  and  Jackson,  it  is  possible  to  answer  the 
stability  question  by  only  evaluating  certain  integrals.  As  shown  in 
Appendix  V,  the  number  of  integrals  equals  the  number  of  maximum 
and  minimum  of  the  distribution  function.  For  a  zero-order  dis¬ 
tribution  function  with  one  maximum,  the  stability  criterion,  from 
Appendix  V  is 
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where 
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and  ^  ,  p  are  the  values  at  the  maximum.  Since  3*^1  =  o  the 

a  rhe 

integral  is  not  singular. 

For  a  Maxwellian  distribution 
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we  find  from  Appendix  V,  that  A  f  -  fZ  .  This  together  with  eq.  40 
gives  the  stability  criterion. 


31 


For  a  resonance  distribution 


t 
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it  is  possible  to  solve  the  set  of  equations  39  exactly.  The  details  are 
in  Appendix  VI.  The  result  is 


*  h,.cr.\  2-44 

V*  / 

When  1  is  negative  then  i  ->  - 1  .  The  criterion  for  stability  is  simi¬ 
lar  to  that  of  the  Maxwellian  distribution.  The  physical  meaning  of  the 
stability  criterion  is  the  following.  For  a  cold  beam  there  is  instability 
if  the  particles  have  a  negative  effective  mass.  This  means  if  u  >  o 
Suppose  then  that  this  is  true.  Then  the  growth  rate  of  the  instability 
is  given  by  the  second  term  on  the  R.H.S.  of  eq.  44.  The  first  term 
on  the  R.H.S.  gives  the  spread  in  angular  velocity  of  the  components 
of  the  beam  due  to  the  spread  in  pe  .  The  stability  criterion  means 
now  that  if  the  spread  in  angular  velocity  is  greater  than  the  growth  rate, 
there  will  be  stability,  because  the  particles  will  have  mixed  themselves 
during  the  characteristic  growth  time  so  that  any  perturbation  will  have 
been  washed  out.  This  effect,  where  a  finite  temperature  effects  sta¬ 
bility,  occurs  in  many  plasma  physics  problems  including  for  example 
the  two- stream  instability. 
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Eq.  44  also  shows  that  when  there  is  stability,  the  oscillations 

are  damped.  This  effect  has  not  been  previously  noticed.  The  result, 

eq.  44,  without  the  damped  solution  is  similar  to  the  result  of  Kolo- 

(9) 

menskii  and  Lebedev.  This  damping  is  due  to  the  additional  term  in 

the  second  of  eq.  2-39,  which  in  the  case  of  one-dimensional  plasmas 

(18) 


gives  Landau,  damping 
because 


For  a  pulse  function  this  term  is  zero, 
<8> 

“  -  O  ,  and  therefore  N-S  who  used  pulse  functions 


3  f 

did  not  find  this  damping.  We  note  that  also  in  the  case  of  one-dimen¬ 
sional  plasmas,  the  use  of  pulse  functions  leads  to  no  damping  and,  in 
fact  gives  the  fluid  equations  for  longitudinal  oscillations,  which  are 
known  not  to  exhibit  damping. 

Finally  we  use  the  results  of  Appendix  VI,  which  gives  the  sol¬ 
ution  of  the  dispersion  relation  for  a  pulse  function  of  width  A  , 


-  6>„  *  ±  _  i*  -  r  (K’o  c  roy  2-45 

l  y  '"or/  y/ 


(8) 

This  result  is  similar  to  that  obtained  by  N-S  and  Kolomenskii  and 
(9) 

Lebedev.  If  =  y*  ,  then  the  result  is  the  same.  This  occurs 

only  when  the  beam  is  sufficiently  thin,  for  then  g  gets  large  and 

£5  o  .  It  is  also  true  for  large  wave  numbers  for  then  <fg  -*o  , 

3 

as  is  apparent  from  Appendix  III.  Since,  by  its  definition  in  eq.  28, 
may  be  negative  for  low  wave  numbers,  we  see  by  eqs.  40,  44 
and  45  that  even  for  a  cold  beam  there  is  no  N.M.I.  for  weakly  rela¬ 
tivistic  beams.  As  an  example  take  the  case  j  =  3,  ;  g  sr 
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(£<jj  is  obtained  from  eq.  C-13)  where 

J_  _  ±  _  J_ 

V  ‘  •r 

For  y>7  ,  the  three  lowest  wave  numbers  are  stable.  For  1=1  , 
=1.3.  In  general,  S<]2  is  a  function  of  2  only.  Our  result 
differs  from  the  other  authors  because  the  beam  is  free  and  is  not  in  a 
vacuum  tank  which  would  shield  parts  of  the  beam  from  another.  Another 
contribution  is  that  the  validity  of  the  constraint  equation,  requires,  by 
eq.  7,  that  the  R.H.S.  of  eq.  45  be  much  less  than  y  ~h  ©.  .  Since 

the  beam  minor  radius  is  much  less  than  the  major  radius,  eq.  6,  the 
constraint  equation,  implies  that  the  first  term  in  eq.  45  is  much  less 
than  0m  .  Hence  the  second  term  must  also  be  less  than  6m  , 

If  4=0,  the  instability  growth  rate  is 
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and  is  valid  only  if  this  quantity  is  much  less  than  6P  . 

From  the  above  results,  it  is  apparent  that  the  criteria  for  sta¬ 
bility  are  insensitive  to  the  precise  shape  of  the  distribution.  Thus 
the  pulse  function  distribution,  which  is  simplest  for  computations, 
gives  adequately  accurate  results.  We  suspect  that  this  holds  true 
for  many  calculations  in  plasma  physics,  where  the  utility  of  the  pulse 
function  is  insufficiently  appreciate. 

In  systems  of  azimuthal  symmetry,  p0  is  a  constant  of  the 
motion.  Thus  neglecting  any  azimuthal  instabilities,  in  a  device,  for 
example,  like  the  betatron,  4  f0  in  eq.  40,  may  be  calculated  from 
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A  fB  of  the  particles  at  the  time  of  injection,  as  it  is  constant.  When 
the  external  field  is  zero  f0=  r  •  Since  r  -*•  rQ,  the  spread  in  p^ 
comes  mostly  from  the  spread  of  v/}  .  The  number  of  particles  in 

(•VI.  - 

a  range  dl/J  is  /!/„  =  ^  a  feT"  .  Hence  the  number  of 

particles  in  a  range  d  p.  is  ,  , 

9  h/  =  / - L - 

f  (tt  •  a  kT  hi.r.3- 

and  comparing  with  eq.  42,  =  -yja  h,.  feT  r,a  .  Inserting  into  eq.  40, 

we  obtain 


— .-1  t  (s.O  I/O'*)*  Te 

L,  ■%  .  • 


where  the  electron  temperature  is  expressed  in  electron  volts.  This 

expression  is  also  valid  for  the  ions  if  W  is  replaced  by  vI-JLl  e 3 

»irr. 

Thus  for  a  plasma  at  a  given  temperature,  that  is  non- relativistic,  the 
number  of  ions  or  electrons  that  are  stable,  is  the  same.  The  quantity 
c/.  is  •<  =  — —  -  and  g  is  given  in  Appendix  III.  For  ions  or 

non-relativistic  electrons,  some  typical  parameters  are  c*  -x  1 1 
T  ^  3  <?,  w.  .  Eq.  47  then  gives  ),3.  x/o'6  which  corresponds 

to  1/50  amp  of  relativistic  electrons. 


Sec.  2.  5  -  The  Effect  of  Betatron  Oscillations 
a.  Axial  (z)  Oscillations 

We  shall  next  investigate  the  effect  of  allowing  small  amplitude 
betatron  oscillations  in  the  zero-order  distribution  function.  To  sim¬ 
plify  the  investigation  we  shall  consider  first  the  (z)  oscillations  and 
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then  the  (r)  oscillations.  We  will  show  that  these  oscillations  have  a 
small  effect  and  give  a  negligible  contribution  to  stability. 

This  result  differs  from  the  investigations  at  MURA  as  reported 
by  Nielsen  et  al.^30)  They  define  two  quantities  ^Ep  and  AEZ  , 
which  give  the  spread  in  beam  energy  due  to  a  spread  in  pfl  (  4  Ep)  and 
a  spread  in  betatron  oscillation  amplitude  (  4EZ).  These  quantities 
may  be  obtained  by  expansion  of  the  Hamiltonian,  eq.  A-l,  since  these 


energy  spreads  are  small.  The  quantity  4  f  p  /  and  A  E  -  X? 

Therefore  the  square  root  of  the  L.H.S.  of  eq.  40  is  proportional 
to  AEp.  (  r*  in  eq.  40  to  correspond  with  Nielsen's  equations.) 

Nielsen  et  al.  now  remark  that  the  effect  of  the  betatron  oscillations 
is  to  add  a  4  Ez  term  to  the  4Epterm,  so  that  both  terms 

contribute  equally  to  stability  .  This  result  is  clearly  different  from 
our  result,  given  below,  eq.  57. 


By  neglecting  the  pr  and  r  terms  in  the  Boltzman  equation, 
eq.  1,  following  the  constraint  equation  arguments  given  in  Sec.  2.1, 
that  equation  becomes  after  linearization, 
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This  differs  in  a  few  respects  from  eq.  13.  First,  there  are  two  new 
terms  in  pz  and  i.  Next  t  now  contains  a  factor  describing  particles 
of  small,  but  finite  amplitude  oscillations.  In  addition  6  is  not 
given  by  eq.  16,  but  now  contains  an  additional  term  proportional  to  z2, 
which  results  from  retaining  the  z2  term  from  A0  in  eq.  D-l.  Thus: 
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We  shall  therefore  be  examining  whether  this  additional  variation  in  B 
will  add  to  the  variation  in  Q  due  to  pff  and  hence  add  an  additional 
term  for  stability  in  eq.  39.  We  also  need  the  value  of  the  other  co¬ 
efficients  in  eq.  48.  We  have  from  eq.  A-3b,  for  the  non- relativistic 
case 

foz~e  2-50 
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Keeping  only  first  order  terms  in  the  coordinates  we  obtain  pz 
from  eq.  A-3c, 


h  ~  —  ('  h——°  1 19  = 
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From  eq.  A-2a 
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We  shall  also  assume  that 
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This  will  give  the  same  stability  criteria  as  an  initial  value 
problem.  Eq.  48  now  becomes 
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+  f-i  '^>n  =  O 
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where  the  a1  s  and  b  are  constants  independent  of  pz,  z,  and  are 
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We  note  that  we  are  forming  two  separate  expansions.  One  is 
of  in  powers  of  p  .  The  other  is  in  powers  of  z,  p  ,  p  .  Thus 
since  we  have  kept  only  the  first  order  terms  in  z  and  pz  for  the  last 
two  terms  in  eq.  54  to  be  consistent  we  must  do  the  same  in  the  equa¬ 
tion  for  .  Thus  since  the  zero-order  Boltzman  equation  is 
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we  have 
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Note  that  ^ 
terms  of  eq. 


-  Cl]  i 
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is  assumed  independent  of  (. 
56  are  zero.  The  solution  o: 


and  t  and  therefore  the  other 
eq.  56  is  evidently: 


H-  h  -f  +aj  jy)  where  A  may  be  a  function  2-57 

of  P  »  • 

We  have  now  given  eq.  44  with  the  z  and  pz  dependence  of  all  the  terms 
given  explicitly.  We  also  know  that  the  solution  of  eq.  54  when  the 
terms  approach  zero  is  given  by  eq.  31  with  R  (2)  =  0.  This  is  suf¬ 
ficient  to  solve  the  equation. 
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We  have  ignored  self-field  terms  in  the  expressions  for  pz  and 

z  that  give  first  order  terms  in  the  perturbation  amplitude  in  eq.  48. 

z  has  the  self-field  term  A'z  whose  neglect  may  be  justified  because  the 

use  of  Y,  obtained  from  the  solution  of  eq.  54  gives  zero  for  A'  .  This 

z 

justifies  the  remark  made  above,  that  Iz  remains  zero  during  the 
N.M.I.  The  dominant  self-field  term  in  pz  is  .° ,  which  has  the 
same  z  dependence  as  the  dominant  external  vector  potential,  and  is 
therefore  ignored.  This  means  that  we  neglect  the  transverse  space 
charge  force. 

The  solution  of  eq.  54  proceeds  in  a  straightforward  way  using 
the  method  of  characteristics.  This  solution  is  then  inserted  into  eq. 

C-2  to  obtain  the  dispersion  relation.  These  details  are  in  Appendix  VIII. 
The  solution  obtained  there  for  ^  has  been  verified  by  insertion  into 
cq.  54  directly.  The  result  is  that  t  in  eq.  39a  is  multiplied  by  a 
factor 

/  -  .00  5  I3  h  J)i 

r  ,v 

This  is  equivalent  to  multiplying  y/  by  this  factor.  Examination  of  eqs. 
40,  44  or  45  shows  that  this  has  a  negligible  effect  on  stability.  Thus 

the  inclusion  of  the  axial  betatron  oscillations  does  improve  stability, 
but  negligibly. 

The  Radial  (r)  Oscillations 

Using  the  above  procedure  for  evaluating  the  z  betatron  oscil¬ 
lations  we  shall  now  evaluate  the  effect  of  the  radial  betatron  oscillations 

on  the  N.M.I.  We  shall,  for  simplicity  ignore  the  z  oscillations  in  our 
treatment. 
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The  linearized  Boltzman  equation  is  then: 

- 1  +  6  ^  +  f  +  /V  ^  -  £>  2-58 

c>  f  3  £  3  fe  d  r  e)  /V 

This  is  just  like  eq.  48,  but  now  we  have  substituted  r  for  z.  We  now 
write  the  coefficients.  6  now  has  first  order  terms  in  r  -  rQ,  so  we 
shall  drop  the  second  order  terms.  By  eq.  A-2b 

Expanding  and  keeping  only  first  order  terms  we  obtain 


9  -  0.  4  fo  -  e.(r-r.) 


2-59 


r  ■ 


P«  is  given  by  eq.  50  and  r.like  eq.  52  is  obtained  from  eq.  A -2a,  giving 

^  •  Again  keeping  only  first  order  terms,  we  obtain  pr  from 

eq.  A-3a 


fr  =  f  * 

r* 


r . 


k  (i-h)(r-r.) 


2-60 


To  simplify  the  calculations,  we  will  define 


r~r«  ~  r  -  r„  - 


^  v-o  (/-  h) 


2-61 


Using  eq.  61  to  simplify  eqs.  59  and  6l,  and  using  also  the  $  ,  t  depend¬ 
ence  of  eq.  53,  the  linearized  Boltzman  equation  takes  the  form 

la-*‘>X)t»  -  a.  <*>,„  3 f.  -  x  3fm 

S  fo  'd  fr 


■*  ai  fr  ~  O 


2-62 
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where  the  a's  and  b  are  constants  independent  of  r,  pr  and  are 

a°  ~  ~  i  (fl  ~  1 *  Ik  fo)  a  j  =  M  ui*  (i  -  h  ) 

rS 


b  -  ~  '<■  l  60 

r o 

a,  ~  ell 


a 


hry 


X  -  r-rtt 


2-63 


Since  eq.  62  is  a  partial  differential  equation  in  r,  pr,  we  may 
consider  pft  =  constant  in  eq.  62  and  hence  dr  =  dx.  k  is  defined  in 
eq.  D-6.  We  note  again  that  we  are  forming  two  expansions.  One  is 
of  tA  in  powers  of  ,  and  the  other  is  of  p.  ,  r.  in  powers  of  r  -  rQ, 
P®  *  Pr*  Thus  since  we  have  kept  only  first  order  terms  in  x,  pr  for 
the  last  two  terms  of  eq.  62,  to  be  consistent  we  must  do  the  same  in 
the  Boltzman  equation  for  .  Thus 


ft  ^ 

t  r  =  o 

3r 

2-64 

aa  *  ^  to 

+  a}  pr  3  =  O 

'dfr 

ax 

and 

t  -  A  y  f  x 5  *  a  j  pr3j  is  a  solution  of  eq.  64,  2-65 

where  Ap  may  only  be  a  function  of  pfl  .  Knowing  the  form  of  given 
in  eq.  65,  and  the  dispersion  relation  obtained  when  the  r  -  rQ,  pr  terms 
approach  zero  (as  given  by  eq.  39),  we  can  now  proceed  to  solve  eq.  62. 

In  writing  eq.  62,  the  first-order,  self-field  terms  /lr  /  /)©  A/  (j)' 
in  r  and  pr  which  contribute  to  the  first  order  Boltzman  equation,  have 
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not  been  included.  This  may  be  justified  in  an  approximate  way. 

Consider  first  r.  By  eq.  A-2a,  this  contains  a  term  A'r  which  is  found 

to  lowest  order  by  integrating  Vj  ,  with  a  factor  pr  .  Looking  at  the 

solution  in  Appendix  IX  found  by  neglecting  A'r  ,  we  find  that  A'r  is 

indeed  small,  for  the  contribution  to  A'  is  obtained  only  from  a  term 
’  r 

odd  in  pr  and  even  in  x.  The  largest  such  term  is  the  third  term  in 

P3 

eq.  1-13  or  14,  and  therefore  AL  has  a  factor  -  ,  which  allows 

r  r  a 

us  to  neglect  it.  This  means  that  the  perturbed  vector  potential  is  small 

because  the  transverse  currents  that  cause  it  are  due  to  particles  which 

have  slow  transverse  velocity  relative  to  the  angular  velocity. 

The  other  neglected  self-field  terms  are  in  pr  .  The  largest  of 

these  is  '  .  The  major  contribution  to  this  term  will  be  a  term 

7>r 

proportional  to  x,  which  is  the  same  spatial  dependence  as  the  external 
field.  Thus  since  the  zero-order  potential  ^  ,  is  assumed  negli¬ 

gible,  we  shall  also  neglect  this  term.  This  means  that  we  are  neglect¬ 
ing  the  transverse  space  charge  forces. 

The  solution  of  eq.  62  is  carried  out  in  Appendix  IX  by  the 
method  of  characteristics.  This  solution  has  been  verified  by  insertion 
directly  into  eq.  62.  The  result  is  that  Vo  in  eq.  39a  is  multiplied  by 
a  factor 

i  -  if  H 

e  /-  <•'  r.x 

l  +  1  Y3Sl 

As  is  more  clearly  seen  from  the  solution  1-22,  this  shows  that  the 
effect  of  the  radial  betatron  oscillations  on  stability  is  also  quite  small. 
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The  denominator  is  unity  if  eq.7  ,  the  condition  for  the  validity  of  the 
constraint  equation,  is  true.  This  formula  thus  suggests  that  as  the 
constraint  equation  is  violated,  the  instability  growth  rate  slows  down, 
presumably  because  energy  is  also  being  put  into  the  excitation  of 
radial  betatron  oscillations. 
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CHAPTER  III  -  TWO  STREAMS 

In  this  chapter  the  dispersion  relation  for  the  N.M.I.  of 

counterstreaming  positive  ions  and  relativistic  electrons  is  derived. 

The  same  formalism  is  used  to  derive  the  dispersion  relation  of  the 

longitudinal  oscillations  of  thin  beams,  which  is  different  from  that 

(5) 

derived  by  Finkelstein  and  Sturrock.  To  show  the  consistency  of  this 

formalism,  the  dispersion  relation  for  infinitely  wide  streams  is  also 

(19) 

found  and  agrees  with  that  found  by  Bludman  et  al. 

Sec.  3.  1  -  The  Negative  Mass  Instability  Equations 

In  order  to  find  the  dispersion  relation  for  the  N.M.I.  for  two 
streams,  it  is  necessary  to  write  the  linearized  Boltzman  equation 
for  each  stream.  Using  the  same  approximations  as  made  above,  in 
Sec.  2.1  and  2.2,  where  the  single  stream  is  discussed,  we  obtain 
again  eqs.  2-13  and  2-16.  The  restriction  on  may,  however,  be 
greatly  relaxed  if  we  assume  a  neutralized  beam,  for  then  there  is 
no  zero  order  transverse  electric  potential.  There  will,  though,  be 
a  first  order  electric  potential  due  to  the  fact  that  under  the  N.M.I. 
there  is  a  transverse  motion  of  unequal  amount  for  each  stream. 
However  this  term,  of  first  order  in  the  perturbation  amplitude  (^) 
if  inserted  into  the  constraint  equation  through  eq.  A-3a,  will  give  a 
second  order  term  in  the  one-dimensional  Boltzman  equation  2-13, 
and  hence  is  neglected. 

The  procedure  is  to  solve  the  linearized  Boltzman  equation  for 
each  stream  and  use  this  perturbed  distribution  function  to  obtain 
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the  perturbed  potentials,  A'  and  (p  These  are  inserted  into  the 
equations  for  p^  for  each  beam.  This  results  in  two  equations  in 
the  two  unknowns  +  and  pfl  -.  The  requirement  that  there  be  a 
solution  implies  that  a  certain  determinant  be  zero,  which  gives  the 
dispersion  relation. 

In  stead  of  solving  the  Maxwell-Boltzman  equations  of  Sec.  2.2 
as  an  initial  value  problem  as  done  in  Sec.  2.  3,  we  will  assume  that 
all  the  perturbed  quantities  vary  as  e‘  lAt  .  This  will  give 
us  the  correct  dispersion  relation  for  the  growing  solutions  and  will 
therefore  be  sufficient  to  find  stability  criteria.  The  Boltzman 
equation,  eq.  2-13  now  gives  after  solving  for  , 

_  •  ?/l  - - =. 

ft'n  =  ^  f*  3-1 

Next  pfl  must  be  evaluated  to  first  order.  By  eq.  2-15  we  have 
for  particles  of  charge  e  and  velocity  \T& 

f0  =  e  34a  .  e  3-2 

c  3  6  dO 

By  eq.  2-25,  we  have  for  ,  neglecting  p^  as  justified  in  Sec.  2.2, 

Me  _  eg  _3_  f fv\  ^ 

39  ~TT  j[c 

-  il  A_  f(}fz.  t. 

To  3  0  J  \  C- 


_Q 

Jr, 


t,\ 


1 _ 

C3- 


fe 


3-3 
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Since  AJ,  is  independent  of  pe  and  is  normalized  to  —  , 


obtain 


Mo  , 


,y>*  j.i 


Also 


Since  the  6  dependence  of  all  the  above  quantities  is  e‘le  we  obtain, 

=  *  I e  ^ -  e_|£-  || t)ndf0\ 

~i]  [llf  tlnJfo  3"6 

r  ~‘?eal  *  3  ^  v;  *'Yr° 

We  have  written  q-  =  /  +  y  <j  f  »/±l^  q  anl*  dropped  one  term  be- 

v/-*  "7T-  ■* 

cause  -1  «  /  .  Also  some  subscripts  have  been  omitted. 

For  pa_  we  obtain, 

fon  =  ‘ill*-  1 1  %}nJfo  ~  ep-  1  j  'k'ncJfrj 

j^1J)clfo  -  njt'nJf0'j  3-7 

=  Li*2i(('K-'k.  *  3  -i  w-  t)dfe 

ru  J  [  <j  o-  c  *  5  crT7"  j  ' 

Inserting  the  values  of  kK  and  H  from  eq.  1,  we  obtain  from  eqs.  6 

and  7  after  writing  £  -  ~  u+  Lr-  I  _  j  ~  <]  l/l a 

v  c*  rj-  J  o rc* 
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hr 

r. 

[f04 

d  fo 

00 

1 

m 

/ 

U-? 

9.  -  £_} 
l  ' 

and 

II 

I 

e a9 
r*  i 

(fo.  dt- 
rS  D 

-  (H  0  fo  + 

3  yPfc) 

3-9 

/' 

e  -  J1 
■  i 

*  -f 

1 

Note  that  f9+-  fiV1  a  con0tant  independent  o£  p9  . 

To  simplify  wo  rewrite  eqs.  8  and  9  and  define 


thus 

fo4  =  f04  D<  ~(h()f0.  D. 
fo-=  -(nt)fo4D*  4  7";  i’o- 

Iv¬ 

or 

fo>(  I'  D<)  *  f>0.(l>i)D-  =  o 

3-12 

r«.  ('«)»*  '  h°^° 

which  are  two  homogeneous  equations  in  two  unknowns.  For  a  non¬ 
trivial  solution  the  determinant  of  the  coefficients  of  }  ft  _ 

must  vanish.  Thus  we  obtain 
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jpetj  =  (/-  0*)  (/-  ly  D-)  -  0-  0,  =  O 


or 


I-  D+  -  ±  D_  + 


-l  v/+*r 

•  rr 


Dj  D_  =  o 


3-13 


In  order  to  examine  this  dieper.ion  relation  more  definitely,  the  terms 
in  eq.  13  must  be  explicitly  evaluated.  By  eq.  10, 

H. 

J  fe 


o.  H 


3-14 


0. 

and  by  eq.  D-2  and  D-6,  to  first  order  in 


Pa 


B,  =  eB  «■  — ! _ 

=  <9«-  (k  >  o) 


3-15 


Thus 


it 

o.  -  in  /_»*■ 

r. 


dfo 


&•-  ~  Jl.  ~  k  f>0 

Using  now  a  pulse  function  for  Vie  4>-(  *  *  a  l  . 

foi4-  we  obtain 


0  fo>A. 

i 


D_  =  A^’g 


3-16 


a7rr«  |fl.  ~  (See  Appendix  VII) 

I.  is  clear  from  eq.  15  that  it.  a.  i,  .  mea,ure  of  the  velocity  spread 
in  the  beam.  It  is  therefore  more  instructive  to  write 


48 


A  6 


Q  r  h .  A  -  _  1  A- 


3  f  Irt .  r  ■L  I  I  -  h 


I 

~p 


—  0<  A  - 

3  t  try.  r 


3-17 


or 


T  A8  *  au; 


3-18 


and  obtain 


ot. 


fj  e?  a  ~  ~  - 

a 7t ru  /m.  <ar/  (A 0)J  ' " y) 

C*  v  <$ U. 


3-19 


(^-)J  '  r‘f)J 

•where  V;  =  f  60*  and  W  -  JL.  -JL 


.  The  appropriate  minus  sub- 

a-rrr  h,.c3 

scripts  should  be  appended  to  all  the  appropriate  quantities.  The  gen¬ 
eralization  to  D+  is  obvious.  Here  f-\  ,  and  we  need  only  append 

the  +  subscripts.  Thus,  eq.  13  now  becomes,  after  dividing  through 
by  \)  = 


I  _ 


ca  <*  ♦  9 

*>,+ 


1_  <*-  5 

+■  rj-  T 


(tv.y  -r0£f  (a*)*  “(^-  -r.fj 

[(/♦*)»  -  JJ  ?3] 


3-20 


v  *»i  - 

+  f  «i  * 


This  equation  is  of  fourth  degree  in  J1  .  Since  this  is  true  even  with- 

Z 

out  the  third  term,  that  term  will  only  modify  the  value  of  the  roots  re¬ 


sulting  from  the  first  two  terms  but  not  change  their  number.  The 
third  term  may  be  neglected  if  V  is  small  enough  as  appears  from 
the  discussion  below.  A  plot  of  this  equation  with  the  neglect  of  the 
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/ ' 

/ 


third  term  on  the  R.H.S.  of  eq.  20  is  drawn  in  Fig.  4,  where 
F  (  r“  T )  denotes  the  two  terms  on  the  R.  H.  S.  of  eq.  20.  The 
solid  line  denotes  the  region  of  validity  of  the  equation,  i.  e.  ,  for 

r„  J1  near  i r¥  or  l/l  .  The  dotted  line  would  appear  if  the  equa¬ 
tion  were  taken  seriously  in  the  rest  of  the  domain  of  r0  £ 

Under  the  assumptions  used  in  deriving  eq.  20,  we  will,  show 
that  it  may  be  written  as  two  separate  equations.  There  are  four 
roots,  of  which  one  pair  is  near  V,  ,  and  the  other  near  vi 
These  may  be  complex.  Since  the  beams  are  in  a  magnetic  field  at 
the  same  position  and  6  =  ll  holds,  we  have  rhKC=h>^.  As 
long  as  the  ions  are  far  from  being  relativistic,  t r_  y>  i /;  wiu  hold. 

Consider  this  case  and  the  pair  of  roots  near  i r  .  Eq.  20  then  be¬ 
comes 


1  -  -  _£_a  <*-  9  + 

V  i r  *  y*  y 

\r  ’ 


-  C  1  Ot  ,  q 

Ki  t  •* 


v*-i\ 

yor  /  3-21 


Two  more  terms  in  this  equation  may  be  neglected.  In  de¬ 
riving  the  constraint  equation  the  magnetic  field  due  to  the  beam  cur¬ 
rent  was  neglected.  This  means  that  *2*/.  Since  rh,.c 
this  implies  also  Kg  £  * ,  ,  and  hence  « r  %  /  and  we  can  neglect 

the  third  term  on  the  R.H.S.  of  eq.  21.  Also  since  ^  C  ,  the 

relation  ^9  «.  /  implies  that  the  first  erm  on  the  R.H.S.  of  eq.  21 
is  much  less  than  J-  and  may  also  be  neglected.  Thus  we  obtain 


(r»  ~  -  u;  )a  = 


^  T7  3-22 

hi  f  J 
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which  is  the  same  equation  as  the  N.M.I.  dispersion  relation  for  the 
proton  steam  alone. 

For  the  pair  of  roots  of  eq.  20  near  VI  ,  the  equation  may 
be  similarly  reduced.  Setting  r,  ^  ^  1C  in  eq.  2C  we  obtain 


_L  =  -±L- 

\)  hi 


r  _£_3  c*4  <?  4-  —  °/'  ^  //  -  9  VI3  -  V'jw.cx+jj  ) 


Since  we  may  again  neglect  the  first  term  on  the  R.H.S. 

of  eq.  23.  The  term  in  parenthesis  may  be  simplified.  If  t'b 
!  _ I _ 

then  ^  2  14-  4-  V  '*».  q  ^  /-  • 

1  7»  hT.  h'+ 

Recalling  now  that  o i,  =  _L  -  /  we  may  write 

I-  3  ^  I  .  u:3/.,  _  «/  . 

3  <rc>  ^  f~  / 

so  that  eq.  23  becomes 


1 

y 


_£_a  of.  <j 
* 


i  -  q/ _ i  .  a\ 

*">  (>-.»  I 


3-24 


(^)‘  -(v.  -  r.£)J 


z 


This  equation  differs  from  the  N.M.I.  equation  for  the  relativistic 
electrons  alone  by  an  additional  term  in  the  bracket,  which  contributes 

to  stability  if  h>±  .  If  the  bracket  term  is  positive  then  there  is 
stability  even  if  A  \r_  -  0. 

Eqs.  22  and  24  must  be  consistent  with  the  inequalities  men¬ 
tioned  earlier  which  are  required  for  the  validity  of  these  equations. 
The  validity  of  the  constraint  equation  requires  that  I  t/I  -  r  ~  Si  I 


1 1  -h 


L 

^  by  ecl'  2-7-  By  the  definition  of  *41/^  in  eqs.  17  and  18, 
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and  the  constraint  eq.  2-6,  because  the  beam  is  thin,  one  can  show 
easily  that  also  AVj  «  .  These  conditions  are  satisfied  for  the 

negative  stream  if  /  and  V  9±L-  /  .  For  the  positive 

stream  using  ftn.c  =.  hn+\/+  ,  the  restriction  is  stronger ^  yfn. 

is  required.  * 

Summarizing,  we  find  that  the  N.M.I.  dispersion  relation  for 
two  streams,  is  given  by  eqs.  22  and  24  when 

A l/;  «  1/4  ^  Vq  «  jr  m.  f  y  In.  1/1  =  ^  y 

'  ~  5 

are  satisfied  and  the  electrons  are  relativistic  and  the  protons  are 
not.  Eq.  22  is  the  same  as  the  N.M.I.  equation  obtained  if  only  the 
protons  were  present,  while  the  single  stream  N.M.I.  equation  for 
the  electrons  is  modified  by  the  presence  of  the  protons.  If  the  elect¬ 
rons  are  also  non-relativistic  then  eq.  24  reduces  to  eq.  22  where  the 
+  subscripts  are  replaced  by  -  and  the  validity  of  the  equations  is 

given  by  V±  (  V]  i<  h.-  **  ,  and  h,.  1/;  . 

h,t  c* 

Sec.  3.2  -  Longitudinal  Oscillations  in  Thin  Beams 
a)  Restricted  equation 

The  dispersion  relation  for  longitudinal  oscillations  of  thin 
beams  will  be  derived  below  as  a  limiting  case  of  the  two-stream 
N.M.I.  equations.  Since  the  p  terms  in  eqs.  6  and  7  were  neglected, 
this  dispersion  relation  is  valid  only  for  cool  beams,  and  for  values 

of  r.  near  or  XT  .  This  defect,  equivalent  to  a  limitation 

on  V  ,  will  be  removed  later  in  another  derivation,  where  the  pfl 
terms  are  retained. 


Longitudinal  oscillations  describe  a  motion  in  which  the  part¬ 
icles  move  in  a  line  along  the  average  beam  velocity  without  any 
transverse  motion.  In  the  N.M.I.  situation  the  transverse  motion 
is  described  by  the  constraint  equation, Sr  =  r -r0  =  —  A  an(j 

hence  if  — —  o  t  J  r  O  and  there  is  no  transverse  motion.  This 

means  that  an  infinitely  strong  focusing  field  inhibits  the  transverse 
motion.  Therefore,  one  need  only  let  -=»  O  in  eq.  20  to  obtain 
the  dispersion  relation  for  longitudinal  oscillations.  Hence  there 


exists  only  one  dispersion  relation  for  beams  in  a  magnetic  field  and 

it  depends  on  the  value  of  l-h  .  Since,  in  general  i-h  is  small 

and  finite,  the  examination  of  the  instabilities  of  contra- streaming 
(5),  (7)  B 

particles  u9ing  the  thin  beam  longitudinal  dispersion  relation  is 

equivalent  to  assuming  777  =  O  ,  which  is  never  legitimate.  It  is, 

however,  of  interest  to  derive  this  relation  because  first,  it  will  be 

valid  for  accelerators  with  strong  focusing,  of  such  strength  that  it 

may  be  considered  infinite,  and  secondly,  as  will  be  shown  later,  this 

relation  is  valid  for  thin  linear  beams,  and  has  been  derived  else- 
(5) 

where  enabling  the  consistency  of  our  treatment  to  be  checked. 

Let  now  77  -)  o  in  eq.  20.  Since  only  oW  (defined  in  eq.  17) 
depends  on  l-h  ,  one  obtains  the  result  that 


c*  _  -=> 


when 


O 


3-26 


For  this  case  too,  we  will  show  that  eq.  20  reduces  to  eqs.  22  and  24. 
Eq.  20  is  valid  when  the  inequalities 
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t/;  -  r„  =0: 

*+  c 

\ 

=  c  , 

V-  -  To  — 

oC.  C 

i 

J 

2 

hold,  because  then  by  the  discussion  at  the  end  of  Sec.  2.  2,  the  p# 
terms  under  the  integral  sign  of  eq.  3,  the  expression  for  A'tf  ,  may 
be  neglected.  The  neglect  of  the  ptf  terms  outside  the  integral  sign, 
which  arise  from  in  eq.  2,  implies  that  t  when  the  ions  are 

non- relativistic  and  4i/l  a  S-  ,  as  also  explained  at  the  end  of 
Sec.  2.2.  Applied  to  eq.  22,  these  inequalities  require  that  yj m.  u  j 


When  applied  to  eq.  24,  the  inequality 

I  .  „  ..  M 


vg  ir 


f  a  1 /  SHjl  a 


«  I 


results.  We  note  further  that  to  obtain  eqs.  22  and  24  from  eqs.  21 

and  23.  requires  that  /  ,  besides  tc  /  .We  shall 

to* 

also  assume  that  »in+  and  therefore  J.  ^  /  -  \>q  ,  When 

r  '  7  TTT 

all  these  assumptions  hold,  eqs.  22  and  24  result  where  a  +  are 
given  by  eq.  26. 

In  deriving  the  constraint  equation,  it  was  required  that 

I17!  “  |^  — j—  y*'  Since  here  o  ,  this  inequality  is  always 

true  and  poses  no  restriction  on  the  parameters. 


The  longitudinal  dispersion  relation  for  thin  beams  may  be 
put  into  a  more  instructive  form  by  combining  eqs.  22  and  24  (with 
the  substitutions  eq.  26)  to  give 


3-28 
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which  is  valid  when  A  u;  C  (/;  «.  c  ,  aic^S.  ,  vq  «  /  »/«*,.«  I  and 

[  I  >4  *»  '  *.♦ 

^3  ^  ^  ^  ^  l  recorc*  again  the  definitions, 


<1  »  I  f  -  I  >3 

3  > *  h.  ,r.  3in  ,r.  'j 

For  positive  numerators,  the  two  terms  on  the  R.H.S 


i  - 

>a  c>  Y 
of  this  equation 


are  plotted  in  Fig.  5.  The  domain  of  validity  is  indicated  schematically 
by  the  solid  line. 


This  equation  differs  from  the  one  derived  by  Finkelstein  and 

(4) 

Sturrock,  by  the  inclusion  of  temperature  terms,  but  more  important 
by  the  additional  factor  £  ’  g  ,  multiplying  the  second 

term  on  the  R.  H.S.  of  the  equation.  For  small  1  ,  and  large  Y 
the  ii.  term  can  be  large  enough  so  that  this  factor  is  negative. 

This  will  give  an  instability  if  AS.  is  small  enough,  as  is  evident 
from  Fig.  4.  Thus,  in  strong  focusing  machines  below  transition, 
where  the  above  approximations  apply,  the  lowest  wave  numbers  modes 
will  grow.  This  is  a  new  instability. 

For  a  linear  geometry  which  F-S  consider,  =  o  as  is  evi¬ 
dent  from  Appendix  III.  Thus  when  J.j  =o  ,  and  V  is  very  small, 
this  additional  factor  multiplying  the  second  term  on  the  R.H.S.  of 

eq.  28  becomes  .  As  yj  _>/,  the  factor  approaches  2,  al- 

h\  ♦ 

though  when  =  /  ,  the  equation  breaks  down.  Thus  for  very 

small  V  this  difference  means  that  the  stable  longitudinal  oscil¬ 
lations  near  VL  are  given  by  eq.  28  as 
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The  F-S  equation  has  instead  the  factor  instead  of  .  Also 

their  <J  =  which  differs  slightly  from  our  g  given  by 

eqs.  C-16  and  C-13  because  the  geometries  differ.  We  cannot  com¬ 
pare  stability  criteria  with  F-S  because  eq.  28  is  only  valid  for 
small  V  . 

When  eq.  28  is  applied  to  the  longitudinal  oscillations  of  thin 
linear  streams  that  are  very  long,  it  is  necessary  also  to  replace 
^  by  —•  ,  where  any  beam  distrubance  is  represented  by 

(p-  (fkj^  q  *  This  is  true  from  the  following  observations.  If 

- -  =  o  i  <fr  ;  r  -  f.  •  Hence  r  =  r«  ,  a  constant,  in  eqs. 

i  *"  H 

A-2b  and  A-3b,  and  these  equations  correspond  exactly  to  eqs.  A-2c 
and  A- 3c  if  the  substitutions  >  f°  ~  ft  ,  f«  ~  f  ,  S  t=At 

are  made.  Thus  eq.  28  is  valid  either  for  devices  of  circular  sym¬ 
metry  with  infinitely  strong  focusing  or  for  linear  beams  without  mag¬ 
netic  fiel  ds. 
b)  Exact  equation 

Now  we  present  a  derivation  of  the  longitudinal  dispersion 
relation  where  the  p#  terms  of  eq.  6,  7  are  retained.  (We  shall  also 
assume  that  the  beam  is  sufficiently  thin,  or  that  the  geometry  is 
linear,  so  that  g=  g  and  =  o  .)  These  two  equations  are  now, 
with  the  substitution  of  from  eq.  1, 

m  1A 

T0+ 


r  o 
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Since  eqs.  30  and  31  are  linear  in  p#  the  simplest  assumption  for  p#,A 
(the  unknown  quantities)  which  gives  a  solution  for  eqs.  30  and  31,  is 


fo±  *  a  t  ■*-  fe 

If  we  substitute  this  into  the  equations,  each  equation  is  still 
linear  in  pe  .  Let  all  the  terms  in  eq.  30  and  31  be  put  on  the  left  side 
of  the  equality  so  that  they  equal  zero.  Since  the  coefficient  of  the 
constant  and  the  coefficient  of  P<  must  both  be  zero,  each  equation 
now  gives  two  equations.  There  are  now  four  equations  homogenous 
in  the  four  unknowns  cl,  ,6,  .By  setting  the  determinant  of  the  co¬ 
efficient  equal  to  zero,  we  obtain  a  4  x  4  determinant  which  is  the 
exact  dispersion  relation. 

There  are  two  simple  relations  which  allow  the  terms  to  be 
simplified  after  the  formidable  number  of  terms  of  the  determinant 
are  obtained.  Note  first  that  0£  =  &t  -  kt  f0  .  Then  the  integrals 
which  appear  in  the  determinant  are  of  the  two  forms 


57 


ff  it 

'  dfJf 


/ 


b.~  —  -  k  f 
i  1 


r  J7  dr 

Q  -  J1  -  k  ¥> 
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Let  0,  ■  -y-5  X  to  simplify  the  algebra.  Then,  since  f$Jr  - 
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r  Xj 
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A  -  h  f 

!r*  i  -  r)l±  , 


k  If 


3-32 


A  '  It  f 


-- 

h  J  \-k? 

gives  the  first  desired  relation. 

To  obtain  the  second  relation,  note  first  that 

/f7FJ'  ’  r*C  '  l  *4r 

because  is  zero  at  the  limits  *  m  ,  and  is  normalized  to  AL  . 

JTT 

The  second  desired  relation  is  now 

JLY^r  = 


kf 


v  A  f  <>  ■*  JL 
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Jir 
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k-  a-rr 
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-  k  f 
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/?  •  JTT 
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'bf  cXr 

4  A/ 
h-  air 

/A-  hf 
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where  we  have  also  used  eq.  32.  (\£  5  6,  *  '  — )  • 

Besides  these  two  integral  relations,  the  following  relations 
are  needed  too, 

eY - l - ._JL  ill—  1  -  v  wg  I 

f.  (c  zir h.  r.  (c  t’h.r.f  jtt/  <*.  \  ~JT  5“ p 

and 

gJ3.  1  ■  N  -  V*.  3-34 

r«  (c  >n«  r„  jJ  J7T^,  lilt  **  * 

Till  now  all  the  equations  are  equally  valid  for  the  N.M.I.  For  the 
longitudinal  oscillations  one  must  let  ~  o  in  eqs.  32,  33  and 
34  so  that 

~  ~  I  ~  ~  — i  and  /{ *  =_li_  -  ~  1  h  .  ~_o(~  =  -  I 

Many  terms  cancel  so  that  the  dispersion  relation  has  finally  the 
simple  form, 


If  »«  use  a  pulse  function  for  t  ,  as  given  by  eq.  G-l,  to  evaluate 
the  integrals,  eq.  35  becomes 


.'-m 

—  c3t) 

J  f 

/  ’1 

<  J 

3 

h  - r-  f ; 

I1  -K-r  ( 

V.  ■ 

-ran 

t  / 

T- 

(n-r 

where  *  ,r  -  &  -  .  .. 

3M,r.  >  «■  r>  and  there  is  now  no  restriction 

on  the  y  values  allowed  in  this  equation.  There  are,  however,  re- 

fo  if  eq.  36  is  valid  as  written.  Since  none  of  the 


strictions  on 


59 


protons  are  relativistic  t+7  /  and 


1 

C a 

i  -  ^r*  *  A  \r*J  1 

-/«/■♦♦  \J 
(  c  c  r  / 

r 

implies  that  if  l/}  t  that  ( 

►»«r. 

.  Also  if  all  of  the  elect- 

rons  are  associated  with  one  value  of  ¥ 

,  then 

1 

a  ' 

x.  -  a  * 

i-  y- 

c  * 

I  +  *  f  \3 

(  C.  C  /  1  In.  l~o 

j_/|-  M-  ~ 

t*  (  trtn.r. 

'  (-**-  f) 

/J(crhi.  rj  J 

and  it  is  required  that  A  f-  tc  c  ’  These  limitations  apply  also 

*  h,  ■  r.  „ 

to  f6  in  eq.  35,  because  the  validity  of  eq.  D-2  requires  f>e  to  be 
small. 


It  is  easily  verified  that  in  the  domain  where  eqs.  28  and  36 

are  valid  the  equations  give  the  same  result.  The  term  _L  1 

a  *>»♦ 

in  eq.  36.  The  criterion  for  stability  may  now 

be  derived  from  eq.  36.  As  is  evident  from  Fig.  5,  one  need  only 

find  the  minimum  of  the  R.H.S.  of  eq.  36  and  ensure  that  it  is  less 

than  —  .  The  result  is  that  for  »/»it  V  4  t-  — —  _i 

V  J  -  H- 

F-S  give  the  criteria  for  stability,  when  j,  ,  as  vat  jw, 

for  their  equation  which  does  not  contain  the  factors  /  _  (Si  r.W  . 

1  tc  ) 

Thus  eq.  36  gives  an  improvement  in  the  stability  criterion  for 
large  r  . 


hi*  J 


becomes 


I  -  /— ) 
1  u  I 


Sec.  3.3  -  Longitudinal  Oscillations  in  Infinitely  Wide  Beams 

With  the  formalism  developed  above,  it  is  now  a  simple  matter 

to  find  the  dispersion  relation  for  the  case  of  infinitely  wide  beams, 

r  a  7T 

or  beams  with  a  perturbation  wavelength  1±.  or  -  ,  much  less 

l  kf 

than  the  beam  width.  Assume  that  the  beams  travel  in  the  z  direction, 
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and  that  all  quantities  EiiipAi.  have  the  i  dependence 

l  kt  ?  -  i  Jit 


given  by  e 


or 


.  From  Maxwell's  equations 

V-  e  =  wy 

*r1k 


Also 


and 


E  -  -  V(P 

c  2t 


e„  =  -IK?* + 


The  gauge  condition  gives 

V-  4  *  -L  (p  =  O 

c  1 

or 

t  It,  Ak  ~  i  </V  '  ° 

Combining  now  eq.  37,  38  and  39,  we  obtain 


<Pk  '  »»A  -T-, 


'  -  (&)' 

Comparing  this  expression  with  eq.  5,  we  obtain 


3-37 


3-38 


3-39 


3-40 


as  the  charge  per  unit  volume.  Thus  in  eq.  5  one  may  set 

I  I 
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yr^ki  I- 


(ti 
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F rom  eq.  39, 
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h,t 


Si 
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Using  eq.  4,  with  the  p#  terms  included  and  eqs.  32  and  33,  eq.  42 

is  satisfied  if  the  terms  VJ  yin.  q  =  o  and  the  term — — —  in 

“JT  *»■*  /r-  nr 

eq.  33  is  set  equal  to  zero.  The  dispersion  relation  is  now  obtained 
from  eq.  30  and  31,  with  the  proviso  that  the  three  mentioned  terms 
are  zero,  and  g  is  given  by  eq.  41.  Following  the  same  procedure  as 
outlined  in  the  above  section  and  solving  the  4x4  determinant,  eq.  35 
is  again  obtained  for  the  dispersion  relation.  Using  the  value  of  g 
given  in  eq.  41,  the  factor  )  -  ^  r*j  is  now  cancelled  and  we  obtain 


if  ir  e3  N 


-f'K  *wr-  )B"-  £ 


r  J  , 

TW  fm 

Q  -  A  i  f0 

7  - 

L  /•*>".  r/j 


fi  ib  now  normalized  to  1,  on  the  field  ,  since  the  factor 

U 

has  been  factored  out.  Since  —  ~  jr  ~  ”  =  density  of  particles 

in  the  toroidal  geometry,  eq.  43  can  also  apply  to  a  linear  geometry, 

where  n  is  the  density.  Also  one  may  write  ^f0-  Jf.p  f  and  normal- 

/  r°  ' * 
ize  ^  to  1  on  the  field  p^  so  that  each  integral  becomes 

( 

T~  —  I  TIT  “?* 


where  M+  =  h->+  ,  M.~  Making  now  a  final  change  of  variable 

and  letting  ft  =  ,  and  normalizing  ^  to  1  on  the  field  vt  , 

M 

and  letting  1  ,  eq.  43  becomes 

x  x; 

kl  -  ^  iv‘  *  { Tvt 


f  d  Vf. 
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.  /Jl  To  I*  =  J_ 

[~n)  'tfi 


that  appears 


where  *±  C  and  is  valid  when  te  t,<y^cr 

Note  that  the  actual  velocity  of  a  particle  \/  is  given  by  V  ~  V+  *  Vt 

If  the  positive  beam  is  at  zero  temperature  and  moving  with 

zero  velocity,  then  the  first  integral  in  eq.  44  reduces  to  At  and 

AA  .  .  .  .  (18)  J l* 

eq.  44  is  identical  to  a  result  of  Bludman,  et  al.  as  shown  in 

Appendix  X. 

The  reason  for  the  additional  factor 

in  eq.  35  for  thin  beams,  but  not  for  infinite  beams,  as  in  eq.  44,  may 

now  be  noted.  It  arises  because  the  force  between  two  small  elements 

of  a  thin  tube,  far  apart,  moving  with  velocity  r°  ,  is  decreased  by 

the  factor  ,  because  the  electric  field  is  decreased  by  this 

amount.  The  decrease  arises  because  the  force  is  like  that  between 

two  small  charges  and  is  ~±  .  The  longitudinal  electric  field  is 

invariant,  hence  the  field  is  obtained  by  writing  E  -  -A—  ,  where  d 

is  measured  in  the  rest  system  of  the  charges.  Since  in  the  lab  system 

this  distance  is  observed  contracted,  i.e.  ,  Al  =  we  obtain 

<? 

.  The  wavelength  of  the  perturbation  \  =  H.  ,  Cor- 


F 85  a  .* 

^11  2 
responds  to  .  The  forces  and  the  electric  fields  are  the  same 

whether  there  are  actual  moving  charge  clumps  or  a  nearly  stationary 
charge  fluid  where  the  clumps  appear  to  move  due  to  the  phase  velocity 
of  the  disturbance  because  the  charge  density,  only,  appears  in 
Maxwell's  equations.  Thus  the  factor  appears  for  the  thin 

beam.  It  is  also  possible  to  show  from  the  formalism  of  eqs.  30  and 
31  and  the  equations  following,  using  Fa  that  the  longi¬ 
tudinal  field  does  in  fact  have  this  factor  for  the  thin  beam 


case. 
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In  the  infinitely  wide  beam  situation,  the  elements  are  two 
plane  sheets,  the  electric  field  is  E-  3.v<r  ,  and  docs  not  depend 
on  /•  ,  nor  does  the  force,  so  that  the  factor  docs  not 

appear. 

If  only  one  beam  is  present,  then  it  is  possible  to  check  our 
additional  factor  /-  ^ ,  by  making  a  Lorentz  transformation 
from  a  stationary  'beam'  to  a  moving  one.  To  simplify  the  algebra, 
the  beam  is  assumed  cold.  We  shall  assume  the  ion  beam  to  be 
absent.  For  small  currents,  i.  c.  ,  tf.  I  ,  J  ~  Iff2)  p'  '^'1C  rC* 
suiting  equations  derived  from  cq.  36  and  44  are  given  in  Appendix  XI 
as  cqs.  K-8  and  K-l  respectively.  The  details  and  the  Lorentz  trans¬ 
formations  arc  done  in  Appendix  XI  and  the  equations  arc  consistent. 
It  is  suggested  that  the  neglect  of  the  retardation  terms  adds  a  term 
to  g  as  defined  in  cqs.  C-13  and  C-l6,  at  least  in  the 

linear  case. 
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CHAPTER  IV  -  SUMMARY 


The  dispersion  relations  are  summarized  here  (see  p.  123  for 
definitions)  : 


A.  SINGLE  STREAM  N.M.I.  (eq.  2-39) 
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Initial  velocity  distributions  investigated: 
1)  Pulse  function 
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Dispersion  relation  (eq.  G-3) 
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Stability  criteria 
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2)  Resonance  function 


_  A/  A 


37T  7T  fo+AX 
Dispersion  relation  (eq.  F-15) 
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Stability  criteria 


B. 


—  *0  or  4> 
>  a 
J 

3)  Maxwellian  distribution 
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Stability  criteria  (eq.  E-17,  E-18) 
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4)  Any  single  hump  distribution  c  fj  F (ft) 

Stability  criteria  (eq.  E-15,  E-l6) 


where 
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TWO-STREAM  N.M.I.  (eq.  3-13) 
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valid  if  Im  II  >  0 
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2)  Pulse  function 


t.  =  t, 

Dispersion  relation  (eq.  3-20) 
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3)  The  above  equation  may  be  reduced  to  the  two  equations 
(eq.  3-22) 
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Stability  criteria  (eq.  3-24)  At/^ 
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Stability  criteria 
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C.  THE  EFFECT  OF  BETATRON  OSCILLATIONS 

The  inclusion  of  the  betatron  oscillations,  non-re lativistically 


shows  that  the  R.H.S,  of  eq.  A  and  hence  V  in  eqs.  2,  3,  4  and  5 


should  be  multiplied  by 
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Dispersion  relation  (eq.  3-36) 
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The  equation  is  also  valid  for  straight  beams,  if  -%> 

r„ 

The  g  factor,  a  logarithmic  term,  is  then  somewhat  changed. 
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E.  TWO  STREAM  LONGITUDINAL  INSTABILITY  (thick  beam,  eq.  3-44) 


K 


4-12 
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where  \f>>  particle  velocity 
l/*= 

Pi  ■ - 

* 

and  the  disturbance  has  the  behavior 


Ik 

e 


Limits  of  Validity 


a  ^  ^  f  r  if 
A  f).  u  hht.c  r 


p  »  _L 

r  M 
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APPENDIX  I 
The  Betatron  Equations 


We  here  record  the  relativistic  equations  for  the  motion  of  a  charged 
particle  in  electric  and  magnetic  fields  using  the  canonical  formalism. 
The  Hamiltonian  in  cylindrical  coordinates  is: 


Ab¬ 
using  the  six  Hamiltonian  equations  of  motion,  we  obtain: 

±  (*-**') 


A-2a. 


b. 


4  = 

ir(fl  '? 

■  ** )  / '» 

'1M0  -  H_ 

-  «?</> 
c3 

fr  - 

e_£ 

34r 

f  v>'0  _2_  /-  tSL  r  £. 

-Ao)  +  e  i 

.  3±J 

c 

3  r 

3r  (,  r  c. 

'  c 

3  r 

11 

o 

••Sis.. 

ex 

c 

ae 

-t-  e  ^£a. 

c  3  o 

4  e 

c 

94*- 

h" 

ei:. 

c 

94r 

3 

t  <L^£ 

c  3* 

i 

c 

"3T 

If  the  betatron  field  is  azimuthally  symmetric,  then 

- >  A 

sented  by  ,  where 


c. 

-  e  31 
dr 

A-3a. 

-  e 

ire 

b. 

-  e  dip 
di- 

c. 

it  may  be 

repre- 

which,  for  small  r-ra  ,  t  using  IS- V^A  gives  ZJ +- dJnX'-  Also  equation  4 

r.  K  \r  ) 

satisfies  v<g=  v<^y^)  =  o  and  obviously  also  V-ft-  y-  Vx/l  =  £>  .  Equation  4  also 
gives  the  familiar  2-1  condition  at  r  =  rQ. 
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For  the  self-fields  we  have: 


<p(rtt)  =  f  f(r\t-lrzr‘l)dif  fo)  s  (?{?•,  + ICzL'tysr 

I  lr-r'l  )  c  Ir-r-i 


A-5. 


Consider  now  an  equilibrium  situation.  and  the  self-fields 

/}#and  ^  are  assumed  small  enough  to  be  neglected.  The  equation  of  mo¬ 
tion  for  the  2-  direction  is 

ft  -  e  j/o  A-6. 


Equations  2b  and  3c  become 


A-7. 


and 

ft  =  **’  4 


A-8. 


^  r-r.  and  2-  are  first-order  small  terms.  Thus  to  lowest  order 
-G»rt ,  a  constant. 

h\C 

Inserting  these  results  into  equation  6  gives 


A-9. 


With  the  same  assumptions  the  equilibrium  motion  in  the  r  direction 

becomes,  using  equations  2a  and  3a, 

fr=h'r 


Note  that  V-  "Jjo  ,  £_/)6  plays  the  role  of  a  potential  for  thcr  and  2 
r  c 

motion,  and  that  \f0  and  V  are  of  opposite  sign. 
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APPENDIX  II 
The  Generalized  Potential 


We  show  here  that  the  idea  of  a  potential  well  applies  to  all  magnetic 
fields  with  O  *  >?  <■  /  .  In  particular  we  will  show  that 

£*  -  2.  Ah  =  A  -  ±At  B-l 

r  c.  r  c 

where  An  is  the  vector  potential  due  to  an  axially  symmetric  magnetic 

field  where  O'-  H  <•  /  near  the  orbit  but  otherwise  arbitrary.  Pj  -  f0  - 

and  ft  correspond  to  the  value  of  f0  at  the  equilibrium  orbit  r*r„  .  Thus 

eq.  A-10  describes  the  radial  motion  if  f0  is  replaced  by  ,  and  the  r 

and  ?  motion  may  be  described  for  these  arbitrary  magnetic  fields  by 

the  potential  V  -  ~  fl  tJLAa  . 

r  c 

Proof:  We  have 

B  =  8, *  =  *  i  a,ld 

Th“8  A^Xt'l.rdr 

r  Jo 

=  ±(r°  H^rc/r  *■  f  C  B+rJr 
r  10  '0 

and  in  particular  if  is  given  by  eq.  2 

<A,  t  hflfsT  _ 


B-2 


An* 


a  -  n 


a  -  *> 


B-3 


Since  the  expression  for  in  eq.  2  has  been  used,  which  is  valid  only 
near  r„  Aq  in  eq.  3  is  also  valid  only  near  r„  .  The  equilibrium  orbit 
is  found  from  eq.  A- 3a,  with  fr  *  r  ■*  o  .  Since  in  the  equilibrium 
situation  A?  Al  “  0  ,  and  A*o  t<p*  are  neglected,  the  equilibrium 

orbit  is  obtained  from 

9  ( fo  -  \  =  o 

l  r  c  ' 


o 
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or 


fo  _  <?  Mo  _ 

P  7  J7 


Therefore 


if!  *  -  i  2^' 

r J  dr 


=.&</.  -  £_Bo 


B-4 


r  =  r„ 


using  also  cq.  3. 

It  is  now  possible  to  evaluate  £o  _ 

r  c 

ps =  r»-  a. 

and  the  expansion 

J_  =  _J_  -  h  (r  ~r»)  +  h(hti)  (r  - r„  )a 


The  substitution 


r. 


r 
»  »» 


B-5 


B-6 


r." 


are  needed  together  with  eqs.  3  and  4.  Thus 

fa  .  eAft  =  Pi  *  fa  -_e  f0*  <-  6,  r."  / r''*1 
r  c  r  c  [  r  3-m  \, 

=  h  -  S-krL  ([>->') ±  +  r.h-J  r'-) 
r  c  3 'h  ^  f  / 


fl 


B-7 


-  e  B.r.  f  i  t  (i-h)(r  -roy\ 

r  |  3  r  *  J 

Ah  as  given  by  eq.  3  is  incomplete.  Since  <"£-o  there  is  also  a  field 

component  fir  -  "  .  This  term  is  -  e.  /L  n.  t,  .  When  this  term 

3 1  c  3r,' 

is  added  to  the  R.H.S.  of  eq.  7,  we  obtain,  by  comparison  with  eq.  A-4, 


A  -  £  A„  --  Pi  -  £  /!« 


Q.E.D. 
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APPENDIX  III 


Equations  for  the  Self-fields 


By  eq.  A-5 

(p(r‘o[i'j)  =  ef  ,+  ~  i 

'  )  jr -r'l 

We  shall  show  below  that  when  the  minor  beam  radius  is  much  less  than 
the  perturbation  wavelength,  that 

cp  ,fr  ,f.  .ft)drd\  dfrdfidfr  C-2 

We  write  f  as  a  fourier  transform  as  indicated  in  eq.  2-18.  For  a 
typical  term  the  right  hand  side  of  eq.  I  is 

•  ~ •-  -■) 

dr  dt  cle  • 
dfrdfkdf9  C-3 


r  = 


;  10  - ;j2t  *  «'  £  Ir-r'l 


LtlVri  ,A) 

hr -r'l 

■  ( A  till  "-r1""'1  Je 

'  Ir-r'l 

where  A  is  independent  of  &  .  Consider  now  the  physical  situation.  There 
is  a  thin  beam  of  approximately  circular  cross  section  with  a  charge  den¬ 
sity  which  varies  slightly  over  the  cross  section.  There  is  a  sinusoidal 
variation  of  charge  as  6  goes  through  2 w,  which  gives  a  corresponding 
variation  in  potential.  Since  it  is  that  we  wish,  it  is  possible  to  make 
a  few  approximations.  Since  the  beam  is  thin,  the  potential  varies  little 
across  the  beam.  Thus  the  observation  point  P  may  be  chosen  at  the  cen¬ 
ter  of  the  beam,  i.e.,  r '  =  r0  .  Since  the  potential  is  approxi- 

rrately  constant  over  the  cross  section  we  shall  also  neglect  the  transverse 
motion  implied  by  the  constraint  equation,  eq.  2-6.  Thus 

| r-rl=  |r*  +r,a  -arr0  ooa_(e-e'J  +  ^^ifrrT  j/  -  Ifunty-e1)'* 


C-4 
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where 

k‘  ‘  -Joi—  «  /  -  (v'-Att’  if  «  / 

(r"-| '/r.3  VrJ* 

Since  r*  r,  for  a  thin  beam,  the  0  integral  is  a  function  only  of 

.  The  r,  i-  dependence  of  l/'"'  is  weak.  Thus  the  surfaces 
of  constant  I  are  cylinders  about  the  center  of  the  beam.  Because  of 
this  symmetry,  the  potential  at  the  beam  surface  is  given  by  evaluating 
the  integral  at  (r-rc  }af  =  f>*  and  assuming  that  all  the  charge  within 
the  cylinder  is  located  at  the  line  r  -  r.  a  -  o  .  This  is  similar  to  the 
approximations  made  in  calculating  the  inductance  of  a  thin  coil. 
Therefore 


3 


C-5 


and  eq.  3  may  be  written  as 

,  •/„  .  n  .  .)  /„  //  -  A1  cv  j 

I  -  e  '  -  1 _ '  1  Jo  c-6 

/  .7  ru  [  /  -  /{3  oov  3f o 

We  have  also  set  »  lf  r  »  r„  .  Next  we  make  the  substitution  0-o'  =  w 


which 

gives 

HO 

A  /  *  ( i  /  «A  l 

I-- 

e 

*  f  £  e 

■  n  "“f1 

Since 

con  £  a  f  i  oi.< /in, 

C-7 


I-  e;"Vt  / 


(X  U 


C-8 


cos  (/  ,  the  sinrA  term  gives  zero,  so  that 

77  7  ;  jy  3  'Jr  /<  %  I5'-* 

c<)0  /./  J  '  'I 

7  I  „  j/  '  /.■  J  t«v  J  _<£  j'/3 

Note  now  that  the  major  contribution  to  the  integral  occurs  when  the  de¬ 
nominator  is  zero.  Thus  the  contribution  of  the  exponent  is  small.  If 

desired,  this  term  may  be  evaluated  by  expanding  i-  y\.. 

dt 


This 
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will  give  a  negative  imaginary  contribution  to  Jl  ,  caused  by  radiation 
damping,  which  is  however  smaller  than  the  growth  rate  for  the  unstable 
case,  as  is  seen  by  evaluating  the  first  term  in  the  expansion. for  the 
exponent  and  using  eq.  G-3. 

The  remaining  integral  is  now  evaluated  by  elliptic  functions.  We 
have  for  l-l  , 

-  /: 


'-IT 


Co-i  d  </. 

(f-  k'c. 

C  <n.  ci< 


*  >!'  — 

>0  c**  *  y » 

-  U  (J  if  *  -^  ’  </>]  d  </>  =  if  f*  (l-  jcf 


C-9 


V5 


where  we  have  written  u  -  tt  ^  q 

■S 


Thus 


5  cot*4"*  d ^  /'T  (l  -  J  ft' d  <p 

*"  %  K  (i  -  k* y)v* 


C-10 


This  integral  may  be  evaluated  easily,  for  o,  /  using  elliptic  functions, 


(25 


(25) 


For  h  --  l,  3 ,  one  makes  the  substitution  ,  and  using  pp.  181-2, 

these  may  also  be  evaluated  in  term  of  the  elliptic  fct.  K,  E.  Since  k.*l , 
we  set  E,  k  =  1  in  these  evaluations.  K  (k)  has  a  logarithmic  dependence. 


After  some  calculation 

we  obtain 

L,  =  K 

il 

1? 

t 

r 

l,  =  r-i 

Lr  K  -  3V 

3  •  sr 

U  --  K'-f 

i  rv. 

-b 

i 

v 

it 

— J 

The  integral  in  eq.  8, 

with  the  neglect  of  the  exponent,  is 

F,  =  3  f 

co-a.  X  o<.  cf  <<. 

1  i 

p  ~  h.* 

C-ll 


C-12 
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Since  c ,  eq.  12  may  be  written  as  the  sum  of  inte¬ 
grals  of  the  form  of  eq.  10,  whose  values  are  given  by  eq.  11.  Thus, 
after  doing  the  algebra, 

F.=  V/C  V(K  -  ±L) 

e,  =  *(*-*)  ?■<-*(«- W)  c-> 

"  tO*':)  T  E,  ~  <f(K  -  Ax  2 -l  *.o1j 

Thus  eq.  8  becomes  :  +  .,1) 


I  *  e  >4  f, 

Jr. 


C-13 


C-14 


Thus  through  eq.  3,  the  right  hand  side  of  eq.  1  is 

R.H.S  eq.  1  =  j  cl  Jl  5L  9  >  V7  JrJi  d yr d  f>t  d £ L  *  LJlt  C-15 

where  <jt  -  £l  .  The  left  hand  side  of  eq.  1  may  also  be  fourier  analyzed, 
as  in  eq.  2-18  .  Since  eq.  1  is  true  for  any  1  and  Jl  ,  the  integrands 
must  be  equal.  Hence  eq.  2  is  proved  and 

3t  =  £l  ,  K  «  JU  fr.  c-16 

Next  we  use  the  above  method  to  evaluate  an  integral  of  the  form, 

e(  J  drdedl  C-17 

/  jr-r'l 

By  the  same  procedure,  eq.  17  may  be  written  in  the  form  of  eq.  7. 

Since  om.  •<  =  e *  »  e  in  eq.  7,  instead  of  the  term  e1**  ,  we  have 
o< 

_E _ j  e _ •  Thus  fourier  analyzing  eq.  17,  we  see  that  it  may 


*■*  f 


be  written  as 


AlJt  =  e 


LI*  j  (r<* >fr,  P°,  fi)drd^dfrJ  fe  dft 

l~4  / 


where  now 


q  =  J.  ,  V,  *  Em 
A  - - 
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Note  that  <jt  >  from  the  concavity  of  the  curve  in  Fig.  6  . 

For  completeness,  we  note  that  if  the  parts  of  eq.  2  and 

18  are  integrated  over,  we  obtain 

f,A  -  'if  *”*(*} 

r°  /  C-20 

A,Jl  ’  e  §,  (  J  “  (f.)  i  f. 


APPENDIX  IV 


Derivation  of  Q  -  -  k-f 4 


We  here  expand  $  in  powers  of  fm  ,  where 

and  At  is  given  by  eq.  A-4.  Since  we  assume  that  f,  is  small  we  keep 
only  the  first  term  in  the  expansion.  Thus  we  write 


0  =  9 


9  *-  /l£  I  •  pa 

Id  Pol  f 


and  we  wish  to  evaluate  M  .  Note  that  Q  is  a  function  of  fo  explic¬ 
itly,  and  also  implicitly  through  the  dependence  of  r  on  fo  ,  because 

r=*^»)  through  the  constraint  equation  r-r„  =  -  A. 

/-*  VlJi  *  A,&° 

from  the  Hamiltonian  Equation  A-l 

h(h'oCf  D-3 

which  gives  **  J  .  Thus 

i£* 

Using  eq.  3, 


li2-  =  _!_  (h  -  1  e  4  ]  3r  +  i\ 

i  1  (  fo  .  e.  a  \ 

^~r[T  c  8 )  . 
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The  crossed  out  terms  are  zero  by  the  constraint  equation.  Setting 
now  fm.0  .which  gives  r-r0  .  and  using  the  constraint  equation. 
jr-itJi.  '  We  °btain  fr°m  eqs*  4  and  5*  (setting  ra?r  ) 

111]  ~  -_L  ■  _  Hf  -  (si  VhjM  f  J_ 

h'’r  \  c  /  (I-*)  eA.j  (  c  J  ^rJ 

D-6 

=  -  1  -  1  *  i 

ra  |  J  c 4  i>i  r  )  TTF1 

=  T~^(~T^T  +  ~p)  ~  "  k  ,  (k>°) 

where  we  have  neglected  the  slight  r  dependence  of  »a  and  also  written 

~  e  A*  ,  Thus  eq.  2  is 

h,c 

e  =  e.  -  k  fB 


Q.E.D. 
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APPENDIX  V 
Nyguist  Diagram 

(22) 

The  Nyquist  diagram  technique  as  used  by  Penrose  is  applied 
here  to  the  N.M.I.  dispersion  relations,  eq.  2-39,  to  obtain  generalized 
criteria  for  stability. 

To  simplify  the  algebra  we  make  the  following  substitutions  in  eq.  2-39. 
We  write 

and  V.  =  JL  -L.  E-l 

O.TT  K 


and  obtain 


where 


/  =  f  _ Tu-  dv 

i r  4  A. 

Z 


K  - 


//  eag 
TrrrTKS 


E-2 


E-3 


Note  that  f $ dv  v,  /  .  We  assume  also  that  /i*  >  O 

•  *o  J 

We  will  now  find  the  condition  that  Si  has  no  positive  imaginary  part, 
as  a  positive  imaginary  part  means  that  any  perburbation  grows  expo¬ 
nentially. 

The  quantity  l  may  be  a  positive  or  negative  integer.  Suppose  that 
i  -  -  izi  .  Then  eq.  3  is 


cl  IT 


E-4 
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For  fixed  2  ,  the  right  hand  side  of  the  equation  defines  a  function  of 
Jl  ,  for  Jl  with  positive  imaginary  part.  Let  us  call  this  function 
•  Thus 

*$)  - 

Consider  now  the  imaginary  A  plane,  upon  which  is  drawn  the  curve  C. 
See  fig.  7  .  As  the  dotted  portion  of  the  curve  goes  to  infinity,  it  en¬ 
closes  the  positive  imaginary  plane.  By  eq.  5,  this  curve  may  be  mapped 
into  a  curve  in  the  W  plane.  Since  the  function  is  analytic  in  the 
upper  half  plane  with  no  poles,  the  curve  C  is  mapped  into  a  curve  D-, 
which  is  also  described  counterclockwise  and  encloses  all  positive  im¬ 
aginary  values  of  A  ,  some  of  which  may  be  enclosed  more  than  once 
by  counterclockwise  loops. 

On  the  dotted  portion  of  the  curve  C,  *  o  .  Thus  to  obtain  an 


explicit  representation  of  D-,  we  must  obtain  the  value  of  2.  for  Jl 
with  a  vanishing  imaginary  part.  By,  e.g.  ,  J.  D.  Jackson  ^  \his  is 


The  curve  D-  described  by  eq.  6,  is  sketched  in  Fig.  8  for  a  Maxwellian 


distribution.  The  curve  is  labelled  with  values  of  real  Jl  .  This  curve 

(22) 

is  from  O.  Penrose.  The  interior  of  the  curve  encloses  values  of 


positive  imaginary  A  .  For  distributions  other  than  Maxwellian  the 


curve  is  asymmetrical,  and  may  have  additional  loops,  so  that  positive 
imaginary  values  of  A  may  be  enclosed  more  than  once. 


If  7  in  eq.  3  is  positive,  then  the  right  side  of  eq.  3  defines  a 
function  , 


E-7 


The  curve  C  is  mapped  onto  the  W  plane  by  the  function  .  Again 
the  dotted  portion  of  curve  C  corresponds  to  0  ,  so  that  we  must 
find  the  value  of  for  vanishing  imaginary  part  of  fl  ,  in  order  to 
find  the  explicit  expression  for  D+  on  the  plane  w  .  Thus 

dv  E-8 


The  curve  D+  on  the  w  plane  Is  also  counterclockwise  and  encloses 
•  n  the  values  of  positive  imaginary  R  ,  enclosed  by  the  curve  C  in 
the  J1  plane.  For  a  Maxwellian  distribution,  o,  any  symmetrical 
distribution  of  V  ,  the  curves  D-  and  D+  are  identical.  For  asym¬ 
metrical  distributions  the  curves  arc  related  as  follows:  If  wo  replace 
-A  by  -71  in  eq.  6,  we  still  get  D-  but  traversed  in  the  opposite  di¬ 
rection.  The  function  2.  is  now  like  2,  except  that  the  imaginary 
parts  have  opposite  sign.  Thu.  if  we  now  reflect  the  curve  D-  about 

ihe  imaginary  axis,  we  will  get  Dt,  traversed  in  the  correct  counter- 
clockwise  rotation. 

We  can  now  establish  the  generalised  stability  criterion.  The  dis- 

per sion  relation,  cq.  3  is  7  =  -  /  mLtle  • *  .1  .  , 

4  c  '  •  Thus  if  the  point  -)  lies  inside 

the  curve  D,  then  we  have  instability.  Thus  for  stability 


?.  >  -  / 


E-9 
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where  is  the  left-most  point  on  the  boundary  of  D  on  the  negative 
real  axis,  (see  Fig.  8  ),  because  then  the  point  2  =  -  I  will  occur  out¬ 
side  the  curve  D.  (In  Fig.  9  this  left-most  point  can  be  either  point 
1  or  2  and  this  case  will  be  covered  below.)  At  this  point  is  real 
and  negative,  hence  the  imaginary  term  in  eq.  6  and  8  is  zero.  Thus 


|  -  O  for  .+  ? 


E-10 


Since  the  imaginary  part  of  2  goes  from  positive  to  negative  as  we 
traverse  the  point  in  the  counter-clockwise  direction  on  the  curve 
D,  that  is  increasing  values  of  Jl  ,  the  solution  of  eq.  10  corresponds 
to  a  maximum  in  f. 

Inserting  the  solution  of  eq.  10,  IT.  =  -  A  into  eq.  6  or  8,  now  gives 
for  the  negative  real  part  of  , 

J£ 


3  (  Jil 

/  v  —i 


u 


\r 


E-ll 


Since  eq.  10  is  true,  eq.  11  may  be  transformed  into  an  integral  without 

,  ,  (22) 

the  principal  part,  (see  O.  Penrose  ) 


p 


2L  . 

3  w  di r 


v  -  v* 


(jLl £»  ^ 


E-12 


This  last  expression  is  finite  because  {  r  f «  +  Inserting 

eq.  11  into  9  we  obtain  finally 


/ 


j°_l£  clv 


=  o  fill)  <  o 


E-13 
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as  the  condition  for  stability,  for  a  single  peaked  distribution.  T.he 
left  side  of  the  above  inequality  is  positive,  and  implies  that  if  the 
width  of  the  velocity  distribution  is  large  enough  the  beam  is  stable. 

(if  t?  <  o  in  eq.  3,  then  the  /  in  the  inequality  is  replaced  by  -/ 
and  the  inequality  is  always  true,  so  even  narrow  beams  are  stable). 

This  condition  is  necessary  and  sufficient  only  for  single  peaked 
distributions,  such  as  in  Fig.  8  .  For  other  distributions,  as  in 
Fig. 9  ,  it  is  possible  to  locate  all  the  points  on  the  axis  and  by  de¬ 

termining  whether  they  correspond  to  maximum  or  minimum,  one  can 
determine  the  sense  of  the  curve  D  running  through  them.  This  is 
sufficient  to  determine  whether  or  not  the  region  between  any  two 
points  lies  inside  or  outside  the  curve.  One  simply  draws  any  counter¬ 
clockwise  curve  connecting  the  points  in  any  desired  order,  but  such 
that  the  sense  of  the  curve  is  correct.  We  have  not  found  any  theorem 
proving  this,  but  the  reader  may  easily  draw  any  number  of  figures  to 
convince  himself. 

We  now  rewrite  eq.  13.  Using  the  definition  of  f  ,  making  the 

transformation  to  ^  again,  we  now  redefine  a  new  F'JL*  2ZL  *(»  normal¬ 
ly  N 

ized  to  /  on  the  field  f>6  .  Thus 


I 


k  I  ..ft  ~  F  d 


>  K 


tfr ,ra  (  F*-  F  dp 

+  J(r°  ‘  rf 


v 3  <7i. c* 


> 


E-14 
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We  define 


Then  we  can  write  eq.  14  as 


Y  fkcr)1 


which  gives  the  criterion  for  stability. 

Penrose  has  a  plot  of  •  when  F^--£  =  is 


F  - 


<4  iTfF” 


-JL* 

'  e  •“ 


The  negative  of  the  left  hand  side  of  eq.  15  is  then  the 
this  plot.  To  obtain  the  correct  units,  one  sets  - 
figure,  with  the  result 


or 


-{L.H.S.  Eq.  15)  =  -  JL 

*• 


L  y  ~  I/a  <*• 


E-15 

E-16 

Maxwellian,  and 

E-17 

left -most  point  of 
/  ,  in  Penrose' s 

E-18 


using  the  right  hand  side  of  eq.  15.  This  result  inserted  into  eq.  16 
gives  the  stability  criterion  for  Maxwellian  distributions. 
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APPENDIX  VI 
Resonance  Function 


We  solve  eqs.  2-39  for  the  resonance  distribution  function, 

w  -  J/_  _A_  _ >_ 

°  J7T  ir  fi  t-A* 


F-l 


This  function  is  normalized  to  N,  the  total  number  of  particles  on  the 
field  B t  .  To  simplify  the  algebra  we  use  eqs.  E-l  and  E-3  and  de¬ 
fine  f  =  p  1C  ,  so  that  t  ~  .  Then  we  obtain 

J7 T 


h3A  < _ 

ir  (b.^)3  *(1(6)* 


F-2 


Next  we  write  F’/ef',  k  A  *  a'  so  that  eq.  2-39a  become! 

_a£' 

\  IT 


/■"  di 

I  -  ~  /  ^v'  tii 

4  V+Jl 


F-3 


or 


We  shall  evaluate  the  integral  using  the  residue  theorem.  This  equation 
is  valid  only  if  jl  has  a  positive  imaginary  part.  Thus  the  integral  is 
evaluated  differently  depending  on  whether  ?  is  a  positive  or  negative 
number.  Let  us  first  assume  that  }  >  0  .  Then  the  poles  of  the  integrand 
in  the  upper  IT  plane  occur  only  at  \Tt  =  -  6>.  ♦  l &'  which  is  a  pole  of 
order  2.  We  will  evaluate  the  integral  now  by  closing  the  contour  in  the 
upper  \T  plane.  See  Fig.  10  .  (Closing  the  contour  in  the  lower  half¬ 
plane  is  also  permissible,  but  because  there  is  now  a  pole  due  to  -Jk  ,  the 
evaluation  of  the  residue  is 


more  cumbersome. )  Then 
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/  ( t?„  •*  *>')  c!\r 


f  *»"J 


:  IVt  Res  (upper  half-plane)  F-5 


where 


Res  |  =  i(. 

'v=v-  i 
<  rfv 


0,  ’  ^ 


f'jfG. r  v  iW)3(0a  ru -;*•)' 


Q*V 


\r  =  tr 


Jv  ('/'+£j 

('tb.  ^ 

F-6 

(60  +  <r  t-l  a‘) 

(Q.  +  S  tldf' 

- 

l^  +  A 
(  ? 

♦ ;  A'j1' 

=  =  I*’ 

(r+fflraf  W  *£)*  V/' 

Inserting  this  result  into  eq.  5  and  6  gives 


v  =  v, 


I  - 


a  k  K  a' 


3 tr  a' 


"IT 


(-0.  f«  A'»  Va,j 


or 


fi.  *  i  a'  »  1  k  K  ~  o 

Separating  this  equation  into  real  and  imaginary  parts  and  setting 
^  ,  0-'  y*  +d'  gives  the  two  equations, 

A*  ~  +  k  K  ~  0 

a  i  /)  C  =  o 


F-7 

F-8 


Eq.  8  may  be  satisfied  by  setting  A  or  B  equal  to  aero.  Since  *'>o 
and  we  have  assumed  Sl:>0  ,  B  cannot  equal  aero.  Thus  the  solutions 
to  the  two  equations  are  A  =  o  ,  0  -  ♦  or  =  -  A*  1  V)(T  .  Again 

since  fl\>  o  t  the  solution  is  ^  -  1KK  -  A  '  ,  valid  only  if  V^T>  A'  . 

Thus  the  complete  solution  is 


-4  =  a  +  ; 

i 


(Via  -  6') 


F-9 
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If  1*0  ,  then  the  pole  at  -Jl  is  in  the  upper  half-plane  and  we 

close  the  contour  in  the  lower  half-plane.  The  only  pole  in  this  region 
is  at  ITj  s  -  -  it'  .  The  integral  in  eq.  5  now  equals  the  negative  of 

the  residue  in  the  lower  plane  because  the  contour  is  described  clock¬ 
wise.  One  sees  easily  that  now 
I  _ 


Since  the  integral  in  eq.  5  is  minus  this  and  proceeding  as  above  we 
obtain  again  the  eqs.  7  and  8  where  A'  *  t  B‘-  -1--  4*  •  Again 

B1  cannot  equal  zero  because  Jl;>o  is  assumed  and  hence  A  - O  and 
B '  =  ±  Vte  K  i  or  _flL  =  -  \/K h  .  Again  since  72 ;  >0  , 

JL;  =  valid  only  if  ■{ KlT  >a' 

Combining  this  result  with  that  of  eq.  9,  we  have 


A  =  90 
l 


VK/d  >  4' 


F -11 


and  the  plus  sign  is  valid  if  l  >  0  ,  and  the  minus  sign  for  J  i  O  , 

Next  we  solve  the  dispersion  relation  for  the  damped  solutions. 

There  is  now  an  additional  term  which  must  be  added  to  the  right  side 

of  eq.  4  as  is  evident  from  eq.  2-39b.  Using  the  transformations, 

t  -  — —  h  "f  i  eq.  E-l  and  E-3  and  t{  a  -  6  ,  this  term  is 

air 


!7 rl(Kk)d£‘(-f) 
r 


3Tri(Kk)fa)A‘r~&-+V - 

~rr[(e,  *■  v)3 


F -12 


v 

i 
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We  assume  that  l>o  and  perform  the  integration  as  in  eq.  5.  Since 
y  now  has  a  negative  imaginary  part,  there  is  a  pole  at  v*  =  -A 
in  the  upper  half-plane  besides  the  pole  at  \T,  .  Thus  we  have  a  term 


3 7Ti  Res 


1 

6.  +«/* 

v--.il 


F -13 


additional  to  the  Residue  at  ^  ,  given  in  eq.  6.  Inserting  this  term  into 
eq.  4  we  see  that  it  cancels  exactly  the  term  in  eq.  12.  The  dispersion 


relation  is  now  again  the  same  as  eq.  7  and  8,  with  solution  f{T0  .  B 


cannot  equal  zero  in  eq.  8  because  then  A  is  imaginary  from  e 
from  eq.  7 ,  B=  ±  AtK  or  *  VffE  .  Since  A; 

negative  the  solutions  are 

-4'  -m.  £■■.-*■*{! a 

i  'l 


q.  7.  Thus 
is  assumed 

F -14 


where  the  first  solution  is  valid  for  all  values  of 

Next  we  examine  the  solution  for  l  -  -|Z /.  Since  the  pole  due  to  -A 

l 

is  now  in  the  lower  half-plane  we  close  the  contour  in  the  lower  half¬ 
plane.  The  integral  eq.  5  is  now  equal  to  the  negative  of  the  sum  of  the 
residues.  Thus  we  have  a  term  equal  to  the  negative  of  eq.  14  due  to 
the  pole  at  -A  which  again  cancels  the  added  term  eq.  12,  because  that 
term  has  a  minus  sign  when  }  is  negative  (see  comment  following  Eq.  2-39). 
The  other  Residue  at  \T-  ir3  is  given  by  eq.  10.  Inserting  this  result 
into  eq.  4  now  gives  the  eqs.  7  and  8  where  A'  ,  B'  are  defined  as  above. 
Again,  the  solution  A'  =  0  to  eq.  8,  gives  real  solutions  to  B,  while  the 

solution  B'  =  0  does  not.  From  eq.  7  we  see  that  B'=±{Kk  or  =  4'i^icl 

Jif 

smce  J7;  is  negative  we  see  that  the  negative  sign  for  the  radical  is 
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valid  only  if  i)Kk  c  a'  .  Summarizing  our  results,  the  roots  of  eq.  2-39 
with  a  resonance  function  are 

A  =  G.  +  l  (-  A'  ±  {Kk  )  F_15 

2 

for  positive  l  .  For  negative  2  »  the  c  becomes  . 

When  K  is  negative,  which  will  occur  if  Jj*<  O  in  cq.  3,  eq.  15  is 
still  valid  although  the  derivation  given  above  must  be  modified.  In 
particular,  the  quantities  &t8'=o  nowand  AtA>M0  . 
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APPENDIX  VII 
Pulse  Function 


Here  we  evaluate  the  dispersion  relation  eq.  2-39, 

i.  e.  , 


function  distribution  for  ^ 

4i  = 

Eq.  2-39  now  becomes 


using  a  pulse 

for  -A  i  p.  *  A. 
a  re  x 

otherwise 


or 


W  r° 

_  ea  q 

air  ft 

/ 

A/ 

,  (*•  -i) 

-  k  “fo 

J7T& 

(e. 

k  i/ 

(e. 

r  eAfl| 

b/ 

-  tea 

V  r. 

JlTT  A 

?■ '  ?/  - 

M 

It  A 
x 


(e.  -  - 

IM3  - 

/  U  e3<j  A/  \ 

l  z  / 

(  JL  / 

[  a-Jrr.  J 

Writing  now  k=  — 


i'Kr 


.  and-— ^  =  V 


r  »„.( 


0  '  aO.  r  ±  U. 
I  — 


t  ra 


—  ~  ^9^"  (hi ,cr) 


,  eq.  2  becomes 


G-l 


G-2 


G-3 


Note  now  that  the  second  of  eq.  2-39  also  gives  the  same  result, 
is  is  because  ||-  *  0  only  When  ±Z±  =  ,  A  ,  which  can  occur  only 
^  by  eq.  3.  Thus  eq.  8iis  valid  for)*n  ft  positive  or  negative. 
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APPENDIX  VIII 
Betatron  Z  Oscillations 


In  this  section  the  method  of  characteristics  is  used  to  solve 
eq.  2-54.  The  dispersion  relation  for  the  N.M.I.  ,  with  the  inclusion 
of  the  axial  (z)  betatron  oscillations  is  thereby  obtained. 

Eq.  2-54  is  a  linear  partial  differential  equation  in  the  variables 
i>fi'  The  standard  technique  for  solving  this  equation  is  given  in, 

e.g.,  Cohen's  "Differential  Equations."^  The  solution  is  obtained 
by  solving 

d  ft  =  d  ?  = _ dV‘n 

al  a,  <f'n  -( aa  *  H“l 


The  solution  of  the  equation  obtained  from  the  first  equality  is  + 

T°  ^  thc  °thcr  solutio"  we  use  the  last  equality.  This  contains 
*  •  Using  c,  ,  to  eliminate  ft  ,  and  eq.  2-57,  we  find  that  k*Af(c.) 
and  hence  is  not  a  function  of  g  any  more.  Thus  from  eq.  1  we  must 
now  solve,  after  rearranging: 

d^lA  (cto  *  biJ)  f,n  a.  q>'» 

~r~  +  \ — - -  _  "  -.LL _ Iro 


This  is  a  linear  differential  equation  for  f'*  as  a  function  of  * 
Again  using  standard  techniques  as,  e.g.  ,  in  the  book  by  Cohen,  ^ 
find  for  the  particular  solution: 

_  d 

.  Jf*i  (<.-«>* '*J 


we 


)n  _ 


-  /«- 


{JLhJJlZL.  di 


ifi 3  (C,  -  a, 


H-3 
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Wc  do  not  consider  the  solution  of  the  homogenous  equation  for  which 
p-  o  as  ^  does  not  contribute  to  the  N.M.I.  (Note  that  the  first 
factor  is  independent  of  t  . ) 

(27) 

We  next  integrate  the  exponent.  Since  and  using  Pierce 

#121,  132  wc  obtain: 


£  =  /JL L±-k±l  d  ?• 


a  fljQj 


Wc  can  now  see  that  the  solution,  eq.  3,  gives  the  right  limit  as  the 

i  -ft  terms  approach  zero.  From  eq.  2-55,  .  From  the 

definition  of  C,  ,  ma xc,*a>  Px  .  Thus  ifi  a:  t/>J«  and  the 

'  a,  7  p 

second  term  in  E  is  of  order  .  Similarly  max^cfj^a  =  and  the 

third  term  in  E  is  also  of  order  .  Thus  the  factors  of  b  are  of 

a  r* 

order  JP^  .  To  obtain  the  limit  of  as  f>-*o  ,  wc  may  therefore  let 
k  -*  0  in  eq.  3.  The  integration  is  easily  performed  to  give: 


f 


m 


Q,  f 


in 


d/>o 


i 

a > 


H-5 


which  is  the  solution  we  would  obtain  from  eq.  2-54  setting  all  the  tf  terms 

t 

equal  to  zero.  Thus  we  recover  the  correct  limit  from  eq.  3.  Any  arbi¬ 
trary  constants  which  appear  in  the  evaluation  of  eq.  3,  may  therefore  be 
resolved  by  noting  that  eq.  3  must  give  eq.  5  when  2  o  .  (  is  the 

minor  beam  radius,  while  H  is  a  coordinate  in  the  beam.) 

We  continue  to  evaluate  eq.  3  and  perform  the  integration  of  the  term 
in  the  bracket.  This  may  be  done  easily,  only  if  we  assume  the  exponent 
E  to  be  small,  for  then  we  can  expand  and  write  £  •*=  /*-£'  .  By  eq.  4 

sineb  Lc  I  ,  we  therefore  require  also  a0  «  }  =  .  As  we 
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shall  see  below,  this  implies  that  in  unstable  situations  the  growth  rate 
should  be  much  less  than  the  axial  betatron  oscillation  frequency.  This 
is  well  satisfied  for  current  densities  of  interest.  Thus  writing 


where  is  an  odd  function  of  £  .  This  eq.  6  differs  from  the  brack¬ 
et  term  in  eq.  3  by  a  constant  because  ^  f(l  +  x)d/  and  differs 

by  a  constant,  because  although  two  indefinite  integrals  differ 

by  a  constant.  To  find  this  constant,  we  let  g  o  in  the  bracket  term 


of  eq.  3,  which  is  equivalent  to  letting  i-*o  .  Integrating  gives 

These  other  integrals  are  trivial  and  it  is  seen  that  eq.  6  differs  from  j^J 
by  the  constant  ~  .  We  can  now  evaluate  the  £  dependence  of  ‘f' ,n 

by  expansion  of  the  exponents  in  eq.  3  which  we  denote  -f \  .  Thus 

frW'"  '  fl'1'  l£l^j 
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Note  that  when  2-0,  f“Ui,  ..  that  eq.  5  i.  obtained,  which 
indicate,  that  we  have  correctly  chosen  the  con.tant.  Note  too,  that 

therefore  /3  ns  -»  ,  „  ...... 

°  as  g  o  ,  which  fixes  the  branch  of  this 

function. 

We  have  kept  only  those  powers  of  E  which  give  terms  .  Terms 
which  give  higher  powers  of  a.  have  been  dropped.  Eq.  7  contains 

erms  proportional  to  — ^  •  1-  ,  /  .  The  terms  proportional  to  1 

and^  a.^only,  cancel  out.  If  we  keep  higher  powers  of  E,  the  terms  in 

a'“  ’  etc' '  w°uid  also  cancel  because  ae  f*~>o  f  'iL,  while 

.in'  remains  finite  as  f~tO .  We  have  also  denoted  the  origta  of  each 

term  in  small  scrip,  beneath  it.  Note  that  f  “  is  also  a  function  of* 

through  c.  a  „  8*  .  „,£»  This  c,  was  a  constant  only  in  solving  the 

right  hand  pair  of  equation,  in  eq.  1,  but  become,  a  variable  again,  in 

the  solution  for  f\  This  follow,  from  the  theory  of  solving  eq.  1. 

In  eq.  2,  we  have  used  the  positive  square  root  in  writing  a,Wr 

■  a jiJtHli h  f  ln  eq.  7  is  defined  only  for  positive  *  .  Using 

the  negative  square  root  for  f,  ,  wc  find  that  sin"  also  changes  sign 
(cf.  eq.  4). 

VTo  obtain  the  dispersion  relation  we  must  now  insert  f,„ 

into  eq.  C-2  and  <p  and  ^  represent  perturbation  quantities  if  . 

ce  we  are  considering  only  axial  betatron  oscillations,  but  assume 

the  constraint  equation  linking  r  and  to  be  still  valid,  the  integral 

only  over  cifo  *  J  *  .  Thus  cancelling  ^ 

relation  becomes: 

/  *  earn,  (  IT  )J1  . 

r7~  I  H-8 
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-  / 


is  a  function  of  c,  ,  and  hence  is  even  in  £  .  uiw,'  is 

odd  in  ,  while  ]l  -  ^i/*  is  odd  in  /’g.  .  Thus  the  second  and 

third  terms  in  cq.  7  give  zero  when  integrated.  The  term  in 

2  is  even  in  £  and  apparently  odd  in  ft  .  As  noted  above, 

however,  sin  changes  sign  when  ft  changes  sign  and  hence 

is  even.  Also  because  Vi  is  zero  for  large  p0  all  terms 
in  hn  wlt,lout  the  factor  u0  give  zero  jn  t)ie  ^  integration.  Thus 
the  terms  in  eq.  7  that  contribute  to  the  dispersion  relation  are: 


f'  *  ~ 

'in  a. 


I  i  (iLL'VUi^’f  h  fa ?)3/« i  y* f  _  i>u c, *  Uj^ *j\:' 

\aa>)  3“ia$  3  a  <ij  ifaJiTi 


H-9 


The  smallest  term  in  this  expression  is  -  .If  wc  compute 

,  Q. 

Tin  »n  cq.  7  to  one  more  power  of  E,  we  would  gel  terms  ~a.P3  p*  . 

*  *  *  Cl 

Since  wc  want  the  lowest  order  non-vanishing  term  in  /  ,  it  seems  that 
wc  should  also  keep  a.  f*  terms.  However,  all  these  terms  are  odd 
in  *  and  hence  give  a  zero  contribution  to  the  dispersion  equation. 

We  now  integi ate  eq.  8,  doing  the  J ?  it j‘f  integration  first.  The 
i  ,  dependence  of  l/l  is  contained  in  a  factor  ,  which  is  in¬ 

dependent  of  f0  ,  and  as  shown  in  eq.  2-57  is  a  function  only  of  c,  . 
Suppose  for  simplicity  we  define 

N  '  *i i  t't  ±  “ *  f3 


'« t 


foi 


H-10 


O 


'hi  is  normalized  to  unity  on  the  field  2  ye 
constant.  Thus 


N  is  a  normalization 


A/  f  di  d  }>,  -  / 
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Since  the  area  of  the  ellipse  c(=  qa  y3  is  A  r-  if—  Jy  j  ,  we  find 

•  The  integrations  are  now  performed  most  easily  by  changing 
to  polar  coordinates.  Let 

>Va  -  ajr1  ) 

/  with  the  inverse  J 

-<*L*x  6  -  2  >  transformations)  p7  /-  co<i  © 

faTP~TZrTfyr‘  j 

The  2^  ^  integration  of  cq.  8,  now  becomes  using  eq.  9, 

ft*  f,'A  d^dft 


-*J 


I  f_AL_  fnl£3  +  r3c**-3&  -r^rA^o-  H-ll 

•  r  to-i  &  ■  o  jjl 


V  circle^ 


where 


T  >  -  i*  A'  [aT  O'rSJr  l/We1  1-  -^3(5COIj(5  -  0^1, 

1t  “iuHT,  kJ  4  [r  X -  — a 

0  *7 


c/  © 


_  i1  y 


—  ,pr  jlc  •  (?r^-  ±-f-yr 

i } V a i  6  [J  T  j  "  J.  i 

7  r  Y  771  -_jt  \ 

q  Aa,  / «x  [  3  ix  11/ 

-l\  f*  _L  /JZJ  r  -7 3hPy  (.009?) 

77  H(  (**'  V  /V 


H-12 


and 


I  =  — 

*  a. 


/-  73h  (.oofj 


/-.* 


27) 


T  ,  .  * 

in  the  integration  over  6>  we  have  used  Pierce  #176,201.  We  have 
also  used  the  values  of  b,  aa,  from  eq.  2-55. The  only  difference 
between  eq.  8  and  the  dispersion  relation  of  eq.  2-39a,  is  that  now  the 
number  of  particles  is  multiplied  by  a  factor 


*,=  /  ~  23*  JlL  (•  oos) 

/v 


H-13 
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Thus  the  effect  of  ?  betatron  oscillations  is  to  improve  stability,  as 
expected,  but  only  by  a  very  small  amount.  Note  that  even  though  this 
effect  goes  as  7  ,  it  cannot  be  large  for  small  wavelengths,  because 

the  equation  for  the  potential  (C-2)  and  hence  eq.  8  is  valid  only  if  S£.»f> 
We  have  assumed  above  that  tf,  to  make  E  in  eq.  4  small 

so  that  e  ft  1+  C  .  By  eq.  2-55  this  inequality  implies  that 

/  A  -  lO.  *  ?h  f0j  «  Vh  6a 


//L-m  «-  Ikfoj  <,  I A  -  2  O.j  *■  jl  hfoj 


Eq.  14  is  valid  if 

l-ft  *  lOoi  *  }lh  f>0\  «  Jh  O, 

or 

This  requirement  is  net  for  weak  currents  and  small  temperature  spreads 
in  the  beam.  See,  for  example,  eqs.  2-8  and  G-3. 
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APPENDIX  IX 
Betatron  r  Oscillations 

Here  we  solve  equation2-62  following  the  same  general  pro¬ 
cedure  as  used  in  Appendix  VIII. 

By  standard  techniques  we  must  first  solve 

d  f>r  d  X  d  </<  lJL 

-dj*  /V  at .  (f'Jt  dfa  ~(<ic  +  &<•)</' tn  h 

Tfo 

The  solution  of  the  equation  formed  with  the  first  equality  is 
tfa  XJ  *  « i  fr  =  c,  .  Using  this  to  eliminate  from  the  last 
equality,  we  obtain  a  differential  equation  like  eq.  H-2  but  with  by 
instead  of  L  * *  .  Thus 

dVA  +  (c 

d  X  Vo7(c,  - 

The  solution  of  cq.  2  is: 

-  f  -±LLJk±_tlx. 

V/,Jl  -  a ,(p‘n(?  ^ *'JJ 


/ 


9/'„  C  cfX 


1-3 


CJ, 


1 


in 


dt 

d  fo 

l^C.  -  ti,  <i) 


1-2 


This  is  the  particular  solution.  The  solution  of  the  homogenous 
equation  again  implies  that  "'*o  ,  which  does  not  contribute  to  the 
N.M.I.  ,  but  is  a  trivial  radial  pulsation. 
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This  solution  is  similar  to  eq.  H-3,  but  since  x=  ftp.\  2%. 

'r*>'  9 fo 

is  also  a  function  of  V  and  thus  cannot  be  taken  out  of  the  integral 
signs,  i.e.  by  eqs.  2-6l,  63  and  65 


Qt  y 

3  fo 


I  +  ArAL  .  Za*tl_l± - 

3(  )  (  **>v, .(/-*>) 


Iff.  f(c.)  .  ArMi'.)  (:.*.'**  ) 

K,  • r.(,-h)j 


We  have  reinserted  the  constant  C,  ,  after  the  differentiation. 

(27) 

Using  Pierce  Nos.  121,  129,  we  can  integrate  the  exponent 
in  eq.  3  to  get: 


E  -  c/y  -  -JLt.  Oj  ^  yj)  I_5 

Ha\  (c,-  CIJYJ)  '/«„ «i  \ c,  ajilj 


We  will  now  show  that  eq.  3  gives  the  correct  limit,  eq.  2-39, 
for  the  dispersion  relation  as  the  radial  betatron  oscillation  amplitude 
f  °  ’  From  e(l-  2-63,  aaa,=  ( i-h )  &* .  From  the  definition  of  C,  , 
max  C,  at  j>*.  Thus  and  the  second  term  in 

E  is  of  order  .  The  first  term  remains  finite  as  is  * 

unity,  and  /a./  =  Jl  -  l  'Ik  fa  is  independent  of  x  and  fr 
Thus  the  b  term  in  E  goes  to  zero  and  we  can  get  the  appropriate  limit 
by  setting  b  =  0  in  eq.  3.  The  integrations,  using  eq.  4,  are  straight¬ 
forward  and  may  be  done  in  closed  form.  Note  that  the  first  term  in 
eq.  4  is  independent  of  x,  while  the  second  is  ~  x.  Thus  the  following 
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two  integrals  are  required  for  doing  the  integrations  of  eq.  3: 


h 


f  — * 

?S~  ^  = 


-L 


5s  d  v 


_L 

a. 


■!«—  -dMrC  '  fU  * 
>i«j  Nc  • 


1-6 

f  y  =  a.  X  -  fcsfc,  -  a,  7 J J  #  Uk 

^  a*  *  aa  Qj 

We  have  again  abbreviated  V~  =  Va,(c,-aa^).  Now  that  the  integra¬ 
tions  have  been  performed  over  dt/  the  frt  Ya  dependence  of  C,  , 
is  again  inserted  into  eq.  3.  Thus  T~  s  «i  fr  and  eq.  3  becomes 


%*  -  “■  f' 


in 


3/’o 


4  /|  If  .  -«?a3  foaX-q,^; 
3c  ,  /»,  tra  (/-»,)  a.1  «  a, 


1-7 


That  this  equation  satisfies  eq.  2-62,  with  K«o  ,  i.e.  b  =  0, 
is  easily  verified  by  substitution.  The  second  term  of  eq.  7  will  be 
shown  presently  to  give  zero  when  integrated  over  in  the  dispersion 
relation.  To  obtain  the  dispersion  relation  must  be  inserted 

into  eq.  C-2  written  in  terms  of  the  fourier  components.  We  ignore 
the  ctalfr  integration  or  if  we  wish,  assume  that  t ,n  contains  a 

factor  normalized  to  one  on  the  field  *  ft  .  The  dispersion  relation 
is  now,  after  cancelling  , 

}  ~  11  f  jVa  JfoJrJfr  1-8 

r“  / 
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For  ease  of  integration  we  will  define  f(a>Xa  *  fr)  where  xsr-r* 


as 


1/y  C,  :  a3  (r-r<,y  *  a,  /v  ^  a»/?3 

for 

o 

N  is  a  normalization  constant,  and  just  as  in  the  discussion  following 

eq.  H-10,  N  -  ,  where  /4e  -  TrJ-^p .  Note  that  since 

the  c Irdj’r  integration  must  be  done  first.  The  first  term  gives 


/  --  lil-  f  i^a  .  -i 

r*  /  2/..  a. 


1-9 


which  is  the  same  dispersion  relation  as  eq.  2-39.  Integration  of  the 
second  term  of  T/*  is  facilitated  by  going  over  to  polar  coordinates. 
The  procedure  is  exactly  the  same  as  used  in  the  discussion  following 
eq.  H-10.  We  define 


f>/  -  a3  f!* 

■al*  6  * 


y  «  >  ?,* 

We  can  now  write 


with  the  inverse 
transformations 


X  —  ai.i  6 
u»  yr  -  R  C-e?  O 


3£  _  If  jK  -  _J_  A/  [j/C]  -  £££-^2  ] 

Sc,  “  aft  3c,  <»»  «  3K  30i  ju 


1-10 


as  f  may  be  regarded  as  a  pulse  function  which  is  constant  as  R  goes 
from  0  to  e>  and  is  zero  otherwise.  Also  y  =  r-r0  -—Iji—and  in 

'  s>i r,(i-i,) 

the  r  integration,  p*  is  a  constant.  Thus  evaluating  the  Jacobian 

dr  d  fr  -  [51  R  J  a  d  o 

iai 


105 


and  f  %  a  dr  ti  f\ 


We  note  now  that  is  defined  only  for  positive  values.  To  find  the 

correct  form  of  for  negative  jir  ,  we  must  go  back  to  eq.  3 

and  change  the  sign  of  the  square  root  term.  In  eq.  7  this  results  only 
in  f>r  changing  sign.  Thus  the  sin  and  cos  representations  for  x  and 
fr  are  valid  as  6  goes  through  2tt  and  the  integral  in  eq.  11  gives 
zero  so  that  the  dispersion  relation  is  given  by  eq.  9.  Thus  eq.  7 
gives  the  correct  dispersion  relation,  as  asserted. 

It  is  interesting  to  speculate  on  the  significance  of  the  /  functions 

in  the  second  term  of  'I'm  .  Owing  to  changes  of  f>9  of  the  particles, 

the  radial  position  is  displaced.  This  occurs  for  each  of  the  tubes  of 

particles  centered  on  different  values  of  .  The  distribution 

* 

function  of  a  tube  slightly  shifted  differs  by  a  ^-function  from  the  previ¬ 
ous  distribution.  This  distribution  is  ~  J  coa  &  ,  where  the  ^-function 
is  on  the  surface,  i.e.  .  for  a  tube  of  radius  f>  .  This  shift 

however  does  not  affect  the  density  of  particles  and  hence  does  not  con¬ 
tribute  to  the  dispersion  relation  of  an  effect,  which  occurs  because  of 
9  electric  fields. 

We  shall  now  derive  the  correction  to  the  above  results  due  to 
finite  betatron  oscillation  amplitude.  Since  both  terms  in  the  exponent 


106 


of  eq.  3  are  small  we  shall  now  expand  the  exponential  in  the  same 
manner  as  we  did  for  the  ?,  oscillations  and  then  integrate.  We  write 
e%  /  *  C-  *  £_*  and  integrate  in  cq.  3.  We  shall  first  consider  the 
first  term  in  which  is  independent  of  x.  Thus  using  eq.  5: 


/  ^ E  *  a. 

_  ^  1  J  d  a  y 

y  C , 

f  a - 

JijOj  ic,  J 


/>  X 


a  * 

+  J  ~  u*  b  [*  Ax 


(/■>) 


(fO 


1-12 


U""  'IfFl 

(f-V) 


This  result  differs  from  the  integral  in  the  bracket  of  eq.  3  by 

a  constant.  The  argument  parallels  that  following  eq.  H-6.  If  we  set 

x  =  0  in  eq.  12  we  obtain  zero,  whereas  from  eq.  6  we  have  — 

<». 

Thus  the  value  of  the  constant  is  found. 

Writing  now 


K  ’>■?,»  ff;  <(<.)(,» 

for  the  part  of  ^  due  to  the  first  term  of  ^ Vo 

3  /’c 


,  we  have 


'i  -  e  +  eJ  -  £* 


In 


>(i  •  lf^£L-A 


. 
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*.  i  ]  k  ■'I  ■^v'  1  j. 

» 

1  _  Ax 

t 

/  («j 

o«  J 

(«*  °t)J 

Oj  c>j 

-if. 

•  it 

ya.  31 

q.1  .  o.J  4 

T 

■J’  + 

/3«. 

I3-' 

a  &o  \  b  -J  f  -dl-v-V 

:il  (a,  a ,)a 

5*/-  -/f 

- 

£*•/.  -/■/* 

+  /li  4_ 

*\  Aarr")W  , 

1 

A  (f-> 

/  (*'  ’]  V 

\3I  3 

'/  («> 

3  1  Do 

(a,  a  j 

j*. 

'7j‘ 

+■/-!■  -/)- 

iV  \J  X  4  . 

J/-  « 

i  / ' 

+  ~  1 

(v  /( 

•ri’ 

a, 

'J  c, 

/{' 

In  the  foregoing  calculation  the  last  three  lines  are  the  result  of  going 

to  one  higher  order  of  E  than  in  eq.  H-7.  We  have  calculated  terms  up 

to  i  <*■'(>•  Those  terms  which  are  odd  in  fj.  or  x  will  not  contribute 

a. 

to  the  dispersion  relation,  eq.  8.  Thus  the  only  terms  in  eq.  14  which 
contribute  are: 


f1  -  _i  . 

r  hi 


I  f  fr). 

*  (‘I;  U])3 


1-15 


Next  we  must  calculate  the  contribution  ^  i.  e.  ,  the  part  of 


&  1ft 


if^due  to  the  second  term  of  in  eq.  4  which  is  proportional  to  x. 

ofo 

We  again  evaluate  the  bracket  of  eq.  3  by  expanding  the  exponent  and 
writing  e  C  I  I  f  f  I*  .  Thus  we  need  the  following 


(*(1***3:)**  =  ~  ^ 

r~~~ 


«  J.  cl 


[■I 


J 


expression: 


108 


-  tt*  -f 

— 

(ct. 

-  fZLJ  + 

a. 

(/  Ayf 

-i 

(<t,  \ 

+  <■', 

i 

/.7  ct,  y a 

(lit  o  ,  jVt  O  j 

1 

(  o 

tit  X 
a  t  <*j 


by3 


H 


if 

1-16 


i-  f'irzT 

1  {*,  W,)J 


where  as  usual  1  [q", 

If  wo  now  write 

.  n  *  in  ~  JO)  ( 

hn  s  C('  'V  Dc,‘ 


•^X  ,  { 

a  f  -  a  « > 


B. 

I  /I 


1-17 


then 


/-  r 


(Note  that  cq.  17  reduces 
and  a „  . ) 


'X  (/  l  t: 

■i 
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to  the  second  term  in  eq.  7,  for  small  x,  fr 


Instead  of  writing  out  the  expression  f,"  in  full,  we  will  note 
only  those  terms  which  give  a  non-zero  contribution  to  the  dispersion 
relation,  eq.  8.  F rom  eq.  10  wo  note  that  ~  is  even  in  x  and  pr 
Thus  only  terms  even  in  x  and  /V  will  contribute.  Replacing  c,  by 

f  °3  /'J  in  CC1H-  5  and  16  gives  the  /V  dependence.  Since  all 
the  square  roots  are  positive,  the  equations  are  valid  only  for  positive 
fr  .  Thus  in  eq.  3  we  must  substitute  -\f  for  \f  to  find  the  behavior 
for  the  negative  .  We  find  that  the  N T  terms  change  sign  and  so 
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does  di-* 


—  X  .  Thus  the  sT  and  -aiv»  '  terms  are  odd  in  pr 


f, 


while  is  also  odd  in  x.  Thus  terms  like  the  following  give 

a  zero  contribution  to  the  dispersion  relation,  fr  a i *  '  X,  X,  pr  '*)) 

a  •  - 1  3  ... 

fr  Ai'y  Y,  fr,  X  ''JA"  X  t  fr  X^etc.  Finite  contributions  are  obtained 
from  terms  like  X 3,  j\  ^  (ja-*,  ‘ )‘  and  X  (*i^‘x)fr 

The  expression  E  is  of  order  o0  t  j>  .  We  thus  find  that 

f_X.  Jy  .  . 

—  is  of  order  f 

f  x  £  t 

'  ^ - -  *  1  is  of  order  a,  j> 

J  !  JLL3  ,(x  is  of  order  c»  ,  /,  Q  „  f  \ ' t  f>  J 

All  terms  of  order  a,"  f,  (V-  a  i,j...)in  f*t  must  be  odd  and 
integrate  to  zero,  because  they  must  be  compounded  of  the  factors 

y)  y  or  y)  fr  .  Thus  the  lowest  order  non-zero  term  is  of 

order  />3  .  There  are  also  terms  of  order  af f*  etc. ,  which 

we  will  neglect  because  we  assume  /  .  Thus  to  terms  of 


order 


f 


r  K  -  I  (  *  (>iC  y  x  P  .  /  X 


and  the  non-vanishing  terms  are: 


c 


_  W’ 


o  nt  a 


> « j 


A  N — V 
(<*> «})’ 


a  X  3  *  J  3  j  /V  *  ) 


J  o  , 
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=  --i-  •  <*, 

3  O  J 

Collecting  now  the  terms  in  eqs.  13,  15,  17  and  19  and  inserting  them 
into  the  dispersion  relation/ eq.  8,  we  obtain 


no 


/=  11  \«.2A£  f(c,)  J_(i  «■  ±a 

r«  M  3^  aj\  **  £l  >  « i 


1-20 


+  «,/),  ^  L Ljill  J  fo  J r  J fr 

I  **»  i 3  J 

We  shall  do  the  d  r  d  pf  integration  first.  Thus  is  a  constant  and 
<fr  =  d  x  .  The  first  factor  is 

~  ff  f  (c,)cl*dfr  -  I 

using  the  definition  of  * f '  given  below  eq.  8.  The  next  factor  is 

Ti~  ff  (  ((,)  fv  dv  dpr 


=  PT  _L__  ft  Hi  f?J  /  c«*5 1 

ia>  Vf>'Jo  a,  /a  0  p- 

--  I.1  Z_ 

«3  y 

We  have  evaluated  the  integral  by  changing  to  polar  coordinates  as  de¬ 
fined  in  the  discussion  following  eq.  8.  Finally,  we  integrate  the  last 
factor,  using  polar  coordinates  again,  and  eq.  10.  Thus, 

r,  =//  |f  J.  Jidfr 


3 o  I? UR 


Hi  cto 
<*  i 


1 


~  — f*)  ~  J  *af  \l*Ced0)RdR 

a>  ir/  7,  A  a,  fi  J, 


<*>_dV 

^|CJJ 


Equation  20  now  becomes,  using  Ij,  I^,  and  1^, 


=  -_3 

J 


=  Hh  j  (Hi  .1  (/ 


D  ^(l  ^  e 


«] 


_  1  A,  L 


*2. 


Af  is  normalized  to  N  on  the  field  <9  fe  .  Defining  Ar-±^,  where 
'o  i®  normalized  to  one  on  the  field  ^  ,  and  inserting  the  values 
of  a;and  b  from  eq.  2-63,  we  obtain 


Ill 


I  -  Jl-.11 


V 


avr0 


a  *0'h) 


I- 


afw; 


we  obtain 


3  Vi 


^  fc 
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i  *  1/5 w/  9f°  A'f  -** 

If  we  use  a  pulse  function  for  ^  ,  the  integration  may  be  performed 

as  in  Appendix  VII.  Because  of  the  factor  before  the  integral  in  eq.  21, 
we  obtain: 


i  -  JL-Il 


I 


±sl-i  -3L 


rl  a  $„■*(/- 
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instead  of  eq.  G-3.  Thus  the  betatron  oscillations  have  two  effects. 

One  is  to  decrease  the  effective  number  of  particles,  because  of  the 
/  term.  The  other  shows  that  the  maximum  rise  time  of  the  insta¬ 
bility  is  of  the  order  of  the  radial  oscillation  frequency.  In  the  case  of 
the  z  oscillations,  the  effective  number  of  particles  was  decreased  by 

a  term  ~  ^  .  Note  too  that  IV  Pu  |  in  order  that  eq.  C-2  holds. 

r. 

Thus  this  factor  is  always  small  and  the  finite  radial  oscillations  have 
negligible  stabilizing  effect. 

We  now  show  how  the  additional  term  in  the  denominator  affects 

the  rise  time.  The  dispersion  relation,  eq.  21,  is  valid  only  if  AS  «8„  , 

because  we  have  assumed  a  thin  beam.  For  stability  (A6>Ja  >  y3c-A 

rv*- 

and  hence  for  small  l  ,  the  second  term  in  the  denominator  is  small 
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It  can  be  large  only  when  there  is  instability.  In  the  limit  when 
yq  C_1  »  we  see  that  the  instability  growth  rate  is 

jl  o  ^  ~  (hh) 

l  * 

These  results  are  all  reasonable.  We  expect  the  betatron  oscil¬ 
lations  to  improve  the  stability.  It  is,  however,  a  surprise  that  the 
effect  is  so  small.  That  eq.  G-3  is  invalid  for  large  growth  rates  and 
becomes  eq.  22  is  also  reasonable  because  the  constraint  equation, 

eq.  2-6,  breaks  down  for  large  growth  rates.  Note,  however,  that  in 

,  a.  .. 

writing  eq.  19  we  assumed  that  higher  powers  of  are  negli¬ 

gible.  This  requires  the  validity  of  eq.  H-16,  with  m  i-h  •  For 
the  unstable  case,  as  A  6  «  6.  ,  this  implies  that  the  second 

term  in  the  denominator  of  eq.  22  is  small.  Thus  our  conclusions  above 
about  the  limiting  growth  rate  which  obtains  when  this  term  is  large  are 
not  accurate.  It  seems  likely,  however,  that  the  inclusion  of  the  higher 
order  terms  will  not  change  the  qualitative  result. 
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APPENDIX  X 
Infinite  Beams 


We  shall  show  here  that  an  equation  obtained  from  Bludman  et 

(!9) 

al.  is  the  same  as  one  obtained  from  eq.  3-44. 

In  eq.  3-44  let  to,  represent  a  pulse  function  of  very  narrow 
width  and  let  V;  =  O  .  Then  the  integration  of  the  first  integral  in 
eq.  3-44  is  done  as  in  Appendix  VII.  If  the  width  of  the  pulse  now  goes 
to  zero,  then  this  integral  is  ix  ,  as  t.t  in  eq.  3-44  is  normalized 
to  one  on  the  field  Vt  .  Thus  eq.  3-44  may  be  written  as 

/  -  *  J±-f  aV  c u  J-i 

u*  )  v  -  y 


which  describes  the  longitudinal  oscillations  of  a  relativistic  electron 
beam,  which  is  not  too  hot,  travelling  through  a  stationary  cold  ion 
background.  We  have  also  made  the  change  of  variables  \T-  V_  *  v,  in 
the  second  integral  of  eq.  3-44. 

09) 

Now  we  shall  obtain  eq.  1  from  Bludman  ct  al.  From 

09) 

eq.  2.19,  p.  750 

u*  represents  the  plasma  frequency  of  a  cold  background  plasma. 
This  equation  is  equally  true  if  the  background  has  no  electrons  but 
only  ions.  Then  ‘v*  —r  ,  and 


(v/  3 


k/.. 


+  X 

i  1 


J-2 
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(19) 

By  eq.  2.  2a, 


we’  jB  *  K_JW 
Jl  -  h  <*„ 


i  i 


J-3 


where  we  have  made  the  changes  o'  t  J 3  o  -*  U  *  y 

assume  no  perpendicular  velocity  components  then  \>±-o 

(19) 


tne  equation  on  the  top  right 


v,  fa,  =  l 


.  If  we 
,  and  from 


J-4 


and 


i.  * 


H 


If* 


a  I  r 


cV 

O'' 


/'iv,  3  i/ 


J-5 


-fat  is  normalized  to  /j  ,  the  particle  density,  on  the  field  d  v~  .  Thus 
inserting  eqs.  4  and  5  into  eq.  3  and  integrating  over  the  perpendicular 
components,  we  obtain  ( fi  =  u  -</;/*,.)  , 

y  it  v  u  3  fo  i 

1 -  -  -  d  \f  -  -  t- 

/5hl.  U/-//tf  ~JT 

where  fa  is  normalized  to  one  on  the  field  dtf  .  Also 


\T  Ly 


-^ut 


K 

do  d I/- 


J-6 


\r 


K 

do-  Jl<s 


U/  “  V* 


J-7 


**  3v/ 


'  v  kt 
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Inserting  this  result  into  eq.  6,  and  then  eq.  6  into  eq.  2  gives  eq.  1 
Q(  £•  D« 


116 


APPENDIX  XI 

Plasma  Oscillations  and  Lorentz  Transformations 


We  will  show  that  the  dispersion  relation  for  longitudinal  oscil¬ 
lations  of  a  single  stream  of  cold  relativistic  electrons,  as  obtained 
from  eqs.  3-44  and  3-36,  may  be  obtained  from  the  equations  for  a 
stationary  'stream'  by  a  Lorentz  transformation. 

First  we  will  consider  the  case  of  an  infinitely  wide  beam  as  it 
is  simpler.  Setting  the  number  of  ions  equal  to  zero  in  eq.  3-44,  gives 
=  o  .  Using  next  a  pulse  function  for  */£_  as  done  in  Appendix  VII, 
and  then  setting  the  temperature  term  equal  to  zero  in  eq.  3-44,  the 
dispersion  relation  for  a  relativistic  beam  becomes, 

Jl-  Kfc  *.)*  K-> 

In  a  stationary  system,  i.e.  ,  moving  with  the  beam,  we  see  the 

(28) 

longitudinal  plasma  oscillations,  given  by 

A  =  i  f-Vo 

where 

f^.)J  -  K-2 

is  the  electron  density  measured  in  the  beam  system  and  h\,  is  the 
mass  measured  in  that  system.  The  wave  disturbance  in  the  plasma  is 
represented  by  a  ^  ^  .  ThuB  (  /«.  t  )  form  a  four-vector.  The 

beam  system  is  unprimed,  while  in  the  lab  system,  moving  with  a  ve¬ 
locity  -V  with  respect  to  the  beam,  the  quantities  are  primed.  Thus 

us  =  1  (us1  -  Vk') 
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or  tj  '  =  Jb'  i  V  k 1 

i- 


because  v «  -  1^. 


K-3 


In  the  beam  system,  ^  is  given  by  of  eq.  2  or 

kJi. 


u  -  J \ 


7  v  h.& 


K-4 


We  now  wish  the  value  of  as  measured  in  the  lab  system.  Since  n 
forms  part  of  a  four-vector  we  have  the  following  transformation  equation, 
where  again  the  primed  quantity  is  in  the  lab  system, 

■ 

The  current  j  is  evidently  zero  in  the  beam  system.  Hence  if  *»-/>.  , 
»)'=>-  n.  .  and  by  eq.  4 


Wj  -  1 


'  V  TT  /,  ’e  *  \'4 


K-5 


/  /►»« 


and  therefore  eq.  3  is 


=  1 /k'  ± 


7  n  h  e 


K-6 


This  is  now  identical  with  eq.  1  if  the  appropriate  correspondences  are 
made,  including  where  h'  is  measured  in  the 

l  »«.  /  (19) 

lab  system  as  is  evident  from  the  article  by  Bludman  et  al. 


Next  we  proceed  to  the  case  of  very  narrow  beams,  i.  e.  ,  tie"  I 
,  (29) 

(b  =  beam  radius).  As  shown  by  Sturrock  this  modifies  the  plasma 
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frequency  given  by  eq.  2.  For  thin  beams 


^h3  A  /  _L 

'  (  1> 


=  I*/, 


V-J  •  ^  3,  *  ^ 


K-7 


(V  -  .  i-  >  7  f  u.  /» 


f  J  foh 


This  is  true  non-relativistically  and  hence  holds  true  for  a  relativistic 
beam  if  Ur  is  observed  in  a  coordinate  system  moving  with  the  beam. 

In  the  lab  system  we  observe  a  frequency  Jl.  ,  which  may  be  derived 
from  eq.  3-36.  We  set  the  number  density  of  ions  equal  to  zero,  or 
equivalently  their  mass  infinite,  in  which  case  the  first  term  on  the 
R.  H.S.  of  eq.  3-36  is  zero.  Setting  the  temperature  term  equal  to  zero, 
now  gives 


A*  V.  k'  i  k-8 

if  iU  «  /  •  Also  we  have  set  jj;  t  ,  and  /,*  '  is  measured  in  the 

>  l 

lab  frame.  Wc  now  wish  to  make  a  transformation  from  the  Us  ,  it 
system  to  the  Uj  ,  k '  system  moving  with  velocity  -V.  Wc  have 

=  Y  (  us  >  Vk) 

Y  (Q  •  Vl<)  K'9 


using  also  eq.  7.  Also 

h'=  y(k 

and 


K-10 


-  *£)  . 


k  =  r  (h 


119 


Substituting  this  value  for  k  into  eq.  9  gives 

Uj  '  =  y  (q  f  v)r  (  k'  - 


or 


tJ7  = 


(<?  -*  V)h' 


(Q  *  V)  h  1 


**  '  *  *1  /  *2*. 
r  c  c  c» 

-  Q  k'  *  Vk'  -  Vh'  ( !+  GV) 

l  r  =*  /  4  V 

1  * 

c?|/ 

c  r 

-  Q(hr>‘, 

1  i  <£JL 
c  * 

)k'  t  l7//' 

=  vu' 

4.  <3  It  ' 

using  ^ 

•••a  • 

The  expression  for  Q 

in  eq.  7  may  be  written  as 

= 

'/Trh**  Lll,J  g, 

=  V-  7T  /V  ca/. 

Jn  r  •  Trt3  J 

'is 


K-ll 


K-12 


=  vc>  3s 

This  v;  is  measured  in  the  beam  system  and  is  proportional  to  the  den¬ 
sity  measured  in  the  beam  system.  Calling  this  K  ,  then  since  by 
eq.  5,  h'=  r  ,  we  have  that  V=yK  .  Since  the  condition  for 

the  validity  of  eq.  8  is  «  /  ,  the  corresponding  condition  in  the 

beam  equation  is  p.  a  ,  .  Hence  by  eq.  12,  <?«  c  ,  and  eq.  11  be¬ 

comes  after  substituting  for  Q  in  the  numerator  and  neglecting  Q  in  the 
denominator, 


Us  '  -  V  k' 


+  i 


t.X 


y 


K-13 
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Identifying  this  U>'  with  J\  in  eq.  8,  we  see  that  the  equations 
are  identical  if 


■M.  -  3 

First  ^ s  must  be  transformed  to  the  lab  system.  The  only  quantity 
in  which  transforms  is  fe.  .  By  eqs.  10  and  13,  neglecting  the  V' 
term  which  is  small 

k  r  fu‘  -  Vk‘ yj  ^  jS 


Thus  in  lab  quantities 


V  <3  ,  -  a  _±_ 
\  h‘  I, 


K-14 


while 


=  2 


r. 

7  A 


i  .  13 


K-15 


by  cq.  C-13. 

These  constants  differ  Bomcwhat  because  one  refers  to  a  circular 
geometry,  and  the  other  to  a  linear  geometry,  and  also  there  is  a  factor 
Y  that  does  not  appear  in  eq.  15.  This  suggests  that  cq.  15  is  in  error, 
apparently  because  the  retardation  terms  were  neglected  in  writing 
eq.  G-12.  Thus  since  eq.  14  is  valid  only  for  small  ,  the  factor 

y 

¥  which  should  appear  in  eq.  15  is  probably 
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DEFINITIONS  AND  SYMBOLS 


g 


g 

k 

l 


mi 

m 

n 

N 


Roman  Letter  Quantity  Symbols 
=  vector  potential  of  the  unperturbed  beam,  A®(  r,  0,  z,  t) 

=  vector  potential  due  to  the  perturbation,  A1  (r,  &,  z,  t) 

=  the  external  betatron  potential,  defined  by  eq.  A-4 
=  magnetic  field  at  the  equilibrium  orbit 
=  velocity  of  light 

=  J  (  Jk  ■>  .U  )  for  small  '/  ,  *1  ,  and  2  *  o  ;  for  2  =  0 

I'llf  lllf 

write  l  =.14.  Its  value  is  ^  5.  For  larger  J  it  is  given  by 
eqs.  C-19,  C-12  and  C-13. 


S6Z 

=  gz  +  -L  .  For  2  =1, 
accuracy  see  eqs.  C-19,  C-12  and  C-13. 

=  8t 

,  where  <*  -  — 


~  =  1.  3,  also  gtf  =  g,  .  For  more 


y  in.  r 

0 

=  Sir  ,  the  magnitude  of  the  wave  vector  when  the  disturbance  is 

r  m  e  <  k,i  -  (  fit 

of  the  form  g 

-  integers,  i  1,  2,  3 . it  gives  the  spatial  dependence  of  a  dis- 

i?o  -l  fit 

turbance  through  <? 

=  rest  mass 

=  rest  mass  of  the  positively  or  negatively  charged  particle 
=  relativistic  mass,  equals  )/m0 
=  field  index,  the  exponent  in  Bz  =  B0^I±j 
=  total  number  of  particles  in  the  beam  of  either  specie 
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Pr  ~  radial  canonical  momentum  (see  Appendix  I) 
pe  =  canonical  angular  momentum 
Pz  =  axial  canonical  momentum 
q  =  electric  charge  density 

r  =  radial  coordinate  in  a  cylindrical  lab  coordinate  system 
rQ  =  equilibrium  orbit  radius 

R  =r0 

t  =  t  ime  measured  in  the  lab  system 

~  ro  0„t  ,  the  average  velocity  of  cither  beam 

z  “  axial  coordinate  in  a  cylindrical  lab  coordinate  system 


b)  Greek  Letter  Quantity  Symbols 


04  =  i-h  ~  »  is  negative  in  the  negative  mass  region 

=  — .Try.  ,  the  relativistic  J'  factor 

— j  ~  ~  S')i  .  gives  the  deviation  from  J_ 

*]  *  .  <  ‘  ']  >x 
the  circular  geometry 

J  9/  =  "  ll}  ,  a  small  positive  quantity  for/7/>  0 


due  to 


A  -  range  of  pfl  values  for  which  a  pulse  distribution  function  gives 
non-zero  values 

Av/i  =  r»  A  r  -rf  (y-p-j'j  f  >  a  measure  of  the  beam  temperature 

which  contributes  to  the  stabilization  of  the  negative  mass  instability 

"X  j-  ’  a  mcasure  fbe  beam  temperature  which  is 

effective  for  pure  longitudinal  oscillations 
A  £  =  width  of  the  betatron  well  (see  Appendix  I) 
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,  for  parabolic  wells  given  by  eq.  A-4,  a  number 


^  _  /i  \  &  R 

e*  =  <1  -  n)  — 

which  is  always  less  than  one 

0  -  angular  coordinate  in  a  cylindrical  lab  coordinate  system 

•  n 

0  =  ~J~ — L_  ,  tho  cyclotron  frequency  at  the  equilibrium  orbit  rf 

t  h,.c. 


V  = 


h/ 


a7 rd  hTT^  *  Budker' 8  parameter,  dimensionless  measure 
of  the  lineal  particle  density  {  V  =  1  gives  17,000  amps  if 
V=  c) 

V.El/ 

1/  =  V 

f  =  minor  radius  of  the  beam,  as  shown  in  Fig.  1 

<P  =  electrostatic  potential  of  the  unperturbed  beam,  (f>‘  ( r,  6  ,  z,  t) 

<p'=  potential  due  to  the  perturbation,  if'  (r,  0  ,  z,  t) 

'fo  =  the  distribution  function  of  the  unperturbed  beam, 

&  ,  z,  pr,  pfl  ,  pz,  t) 

^  =  the  distribution  function  of  the  perturbation, 

^  I  Z,  Pyl  P„  ,  P^>  t) 


c)  Mathematical  Symbols 
^  implies  proportional  to 

**  implies  approximate^  equal  to.  This  notation  follows  the  SUN 
Commission's  recommendations  listed  in  Physics  Tcdav,  15  19 
(1962)  7 — - 

P  implies  principal  value 


d)  Dictionary 


N.M.I.  =  negative  mass  instability 


Transition  Enurgy=  for  strong  focusing  accelerators,  1  -  n  is  replaced 
by  kfl  in  the  expression  for  ,  where  ks  >  1.  The  transition 

onorgy  occurs  for  such  i'  than  «  =  0 


10 

0 

D 

0 

0 

0 

0 

0 

0 

0 


c 

0 
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Fig.  1.  The  physical  model.  Electrons  rotate  in  the  indicated  sense 


Fig.  2.  The  vector  potential  well  due  to  the  external  field  AB  and  the 
self-field  A0.  The  lines  of  the  self-field  B°  are  drawn 
schematically  and  encircle  the  torus,  (see  p  13) 


Fig.  5.  Plot  of  the  two  stream  longitudinal  instability.  F  =  the  two 

terms  on  the  R.H.S.  of  eq.  3-E8.  Note  that  if  the  numerator 
of  the  second  term  in  this  equation  is  negative,  then  the  right 
side  of  the  above  figure  is  inverted  about  the  horizontal  axis, 
(see  p  54) 


B 

Fig.  6.  The  ordinate  is  |i.  .  Since  Itl  the  concavity 

of  the  curve  shows  that  q  >  a  ,  (see  p  79) 
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